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PREFACE 

This  Project  RAND  Memorandum  extends  the  study  reported  in  PJM-2292, 
^ermodynamie  Properties  of  Carbon  Dioxide  to  2lt.,000°K  -  With  Possible  Appli¬ 
cation  to  the  Atmosphere  of  Venus.  Whereas  that  earlier  study  assumed  an 
atmosphere  of  pure  carbon  dioxide,  the  present  study  assumes  an  atmospheric 
composition  of  85  per  cent  carbon  dioxide  and  I5  per  cent  nitrogen  by  volume . 

The  calculated  thermodynamic  properties  of  this  mixture,  as  well  as  the 
normal-shock- wave  properties  at  selected  altitudes,  are  presented  as  an  aid 
to  studies  of  the  aerodynamic  phenomena  of  high-speed  vehicle  entry  into 
the  atmosphere  of  Venus . 


V 


ACKNQWIJEDGMEMTS 

The  author  Is  very  grateful  for  the  many  instructive  discussions  with 
J.  L*  Raymond^  S*  H.  Dole^  and  F.  R.  Gilmnjre.  The  assistance  of  Jeannine 
McGann  and  Carolyn  Huher  with  the  calculations  and  graphs  is  also  greatly 
appreciated. 


vil 


SUMM^Y 

In  this  Memorandum  a  model  for  the  atmosphere  of  Venus  is  developed 

which  may  be  useful  for  design  studies  of  early  atmospheric-entry  vehicles 

and  for  wind-tunnel  simulation  application.  This  atmosphere,  derived  from 

many  questionable  assumptions  and  scant  experimental  data,  consists  of  85 

per  cent  carbon  dioxide  and  I5  per  cent  nitrogen  by  volume. 

The  thermodynamic  properties  of  the  derived  atmospheric  composition 

are  presented  over  the  temperature-pressure  range  of  150°K  to  24,000°K-  and 
-4  2 

10  atm  to  10  atm.  To  further  assist  in  aerodynamic  entry  calculations, 
normal- shock- wave  characteristics  of  such  an  atmosphere  are  also  presented. 
The  fact  that  nitrogen  is  present  has  important  consequences  in  the  thermo¬ 
dynamic  properties  and  electron  concentrations  over  the  full  temperature  and 
pressure  range.  The  effects  are  of  the  order  of  the  percentage  of  nitrogen 
addition  to  pure  carbon  dioxide. 

A  graphical  method  of  obtaining  obtaining  oblique-shock-wave  data  and 
a  method  of  obtaining  electron  concentrations  at  low  temperatures  are  in¬ 
cluded  as  appendixes . 
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I.  IHTRODUCTION 

In  ordsr  that  design  studies  may  be  carried  out  concerMng  space -vehicle 
entiy  into  planetaiy  atmos^eres,  it  is  necessary  that  the  atmospheric  coirrpo- 
sition-pressure-density  relation  with  altitude  for  a  particular  atmosphere  in 
question  be  known  to  a  fair  degree  of  precision.  Since  shock-wave,  deceler¬ 
ation,  and  aerodynamic-heating  calculations  are  all  dependent  upon  this  re¬ 
lation,  it  must  at  least  be  bracketed  within  reasonable  limits.  This  study 
is  concerned  with  estimating  this  relation  for  the  atmosphere  of  Venus. 

Unfortunately,  very  little  is  known  concerning  this  atmosphere.  Through 
theoretical  arguments,  spectroscopic,  photometric,  and  radiometric  measure¬ 
ments,  several  theories  have  been  advanced,  no  two  of  ^ich,  however,  appear 
to  be  the  sarne.^^"®)  An  attempt  is  made  here  to  bring  about  reasonable  agree¬ 
ment  between  two  theories  that  are  approached  from  different  viewpoints.  The 
resulting  atmospheric  model  could  then  be  used  as  a  tentative  mean  ^diich  could 
be  coDi)ared  to  other  theoretical  and  experimental  models  until  more  experi¬ 
mental  data  are  gathered  and  a  better  understanding  of  the  Venus  atmosphere 
develops.  The  establishment  of  this  mean  model  atmosphere  and  its  thermo¬ 
dynamic  and  aerodynamic  characteristics  is  the  objective  of  this  study. 

The  first  theoiy  on  the  Venus  atmosphere  ^diich  is  considered  is  that 
of  Dole.  The  partial  pressures  given  in  Ref.  3  correspond  to  a  per  cent 
by  volume  equal  to  the  reported  per  cent  by  wei^t.  Since  the  molecular 
wei^t  of  Ng  and  CO^  are  markedly  different,  this  cannot  be  the  case.  THro 
methods  for  calculating  atmospheric  composition  have  been  employed.  Using 
the  basic  data  on  pp.  h-6  and  the  theory  of  Dole,  the  probable  range  of 
atmospheric  compositions  has  been  recalculated;  a  probable  average  composition 
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according  "to  "this  "theory  has  also  been  calculated.  The  second  nie"thod  em¬ 
ployed  makes  use  of  some  of  "the  available  experlmen"tal  data,  one  of  Dole’s 
basic  assumptions,  and  one  assumption  concerning  "the  pressure-densi"ty  re¬ 
lation  of  "the  Venus  a"tmosphere  belov  the  cloud  layer.  The  assumptions  in¬ 
volved  are  open  to  debate  and  "the  experimental  data  are  highly  uncertain, 
but  "the  remarkable  result  is  "that  the  t"wo  me"thods  yield  quantitative  esti¬ 
mates  of  the  a"bmospheric  composition  in  close  agreement  "wi-th  each  o"ther. 

It  is  "taken  as  "the  basic  statement  of  this  Memorandum  that  these  methods 
do  yield  a  reasonable  quantitati"ve  estimate  of  "the  atmospheric  composition. 
The  composition  an:i"ved  at  is  85  per  cent  CO^  and  15  per  cent  by  vol"ume, 
and  "the  a"tmosphere  is  "then  modeled^  i.e.,  "the  pres  sure -temperature -altitude 
relation  is  derived.  The  "thermodynamic  properties  of  "the  chosen  composi¬ 
tion  are  calculated  over  the  temperature-pressiire  ranges  of  150°K  to  24, 000°K 
and  10"^  atm  to  10^  atm,  respectively.  Finally,  normal  shock-"V7ave  properties 
are  presented  for  selected  altitudes  in  "the  chosen  a"tmospheric  model. 

The  appendices  include  a  me"thod  of  ob"tainlng  oblique  shock-"wave  properties 
from  normal  shock-"wa"ve  properties  and  a  ine"thod  of  obtaining  electron  concen¬ 
trations  at  low  "temperatures.  The  first  is  useful  when  it  is  impractical  to 
compute  and  plot  oblique  shock  charts  because  of  uncertainty  in  the  charac¬ 
teristics  of  a  particular  a"tmosphere.  The  second  is  useful  for  determining 
the  temperat\ire  and  pressure  influence  on  "the  electromagnetic  en"7lronment, 
for  example,  of  radio  signals  to  or  from  an  entry  vehicle. 
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II.  A  TENTA3?IVE  ATMOSPHERE  FOR  VENUS 

COMDUIONS  AMi  ASSUMPTIONS 

IPvro  methods  will  be  used  to  deduce  the  composition  of  the  Venus  atmos- 
phere;  one  is  that  of  Dole,  vhile  the  other  is  based  on  a  few  selected 
e:q)er±mental  iiieas\ireiasnts.  The  first  approach  will  merely  require  that 
sll^t  numerical  refinements  be  made  in  Dole's  paper;  the  trace  of  argon 
mentioned  In  the  paper  will  be  neglected,  but  the  theory  will  be  assumed 
correct.  The  second  approach  will  require  assumptions  concerning  both  the 
validity  of  some  highly  specmlative  experimental  data  and  the  pressure-ten^r- 
ature -altitude  relation  below  the  Venus  cloud  layer.  More  will  be  said  later 
concerning  the  acceptance  of  these  data  and  the  assumptions  involved  in 
choosing  the  final  model.  However,  before  continuing  with  the  derivation 
It  is  necessary  to  state  that  this  is  not  lirbended  to  be  a  rigoaroxis,  absolute 
derivation  of  the  composition  of  the  Venus  atmosphere.  It  is  solely  intended 
to  provide  some  justification  for  assuming  a  particular  model  atmosphere. 

The  discussion  at  the  end  of  this  section  will  shed  some  li^t  on  the  diffi¬ 
culties  associated  with  the  acceptance  of  this  model;  as  for  acceptance  of 
some  of  the  experi mental  data,  the  reader  may  consult  the  referenced  literature. 

First,  some  general  assumptions  pertinent  to  both  methods  are  neces¬ 
sary.  The  constituents  of  the  atmos^ere  will  be  assumed  to  obey  the  per¬ 
fect  gas  law,  p/p  =  RT/M.  The  gas  mLxtiire  will  be  assumed  to  be  of  con¬ 
stant  composition  throu^out  the  atmosphere;  this  neglects  "settling"  of 
heavier  constituents,  diffusion,  jhotodissociation,  and  ionization.  The 
complete  gas  mixture  is  assxjmed  to  be  in  hydrostatic  equilibrium  in  a 
constant  gravitational,  field,  with  centrifugal  force  due  to  planetary 
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rotation  being  neglected.  Now,  as  pointed  out  in  Ref.  1,  all  of  these  as- 
sunrptions  except  the  perfect  gas  law  may  be  far  from  realizable,  especial¬ 
ly  at  high  altitudes.  However,  they  introduce  far  less  xuicertainty  than 
do  the  asstmiptions  to  follow,  and  they  are  fairly  good  assunrptions  for  the 
altitude  range  of  aerodynamic  interest. 

Following  Dole,  the  amount  of  N^  in  the  Venus  atmosphere  will  be  as¬ 
sumed  equal  to  that  in  Earth's,  but  scaled  down  to  make  allowance  for  the 
difference  in  the  sizes  of  the  planets.  Since  the  atmosiiieric  gases  are 
mainly  generated  at  or  near  the  surface  of  a  planet,  relative  masses  of 
in  the  two  atmospheres  will  be  scaled  in  proportion  to  the  surface  areas  of 
the  two  planets.  Moreover,  the  Venus  atmosphere  is  assumed  to  consist  of 
only  COg  and  N^,  the  trace  constituents  being  neglected.  These  two  assumptions 
are  basically  deduced  by  Dole,  and  they  are  basic  to  both  the  numerical  re¬ 
finements  in  Dole's  paper  and  to  the  second  method. 

The  second  method  will  require,  in  addition  to  the  above,  that  the 
atmosphere  of  Venus  e:cist  in  adiabatic  equilibrium  from  the  surface  to  "the 
top  of  the  cloud  layer."  This  is  another  very  basic  and  important  assumption. 

The  numerical  data  consistent  both  with  viiat  is  known  of  Earth  and  Venus 
and  with  the  above  assumptions  are  presented  below.  For  detailed  information 
on  these  data,  the  reader  is  directed  to  the  references.  The  composition  of 
the  Venus  atmosphere  will  then  be  calculated  from  considerations  of  these 
two  methods  or  approaches  and  the  data  below. 

BASIC  D.ATA 

(9 ) 

A.  Com^sition  of  air  by  volume  for  Earth 
Og  =  20.946  per  cent 
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Ng  =  78.084  per  cent 
+  trace  eleinents 

B.  Molecular  wei^ts  (chemical  scale 

Mj^  =  28.016  g/g-mole 
2 

M^O  =  44.011  g/g-mole 
=  28.967  g/g-mole 

C.  Specific  heats  (450°K)^^^^ 

(C  =  7.017  cal/g-mole  °K 

(C  )„  =  5.031  cal/g-mole  °K 
(Op)co  =  10.280  cal/g-mole 
(C^)co  =  8.294  cal/g-mole  °K 

D.  Mass  of  COg  per  vinlt  area  of  Earth's  surface  (in  the  atmosphere, 
water,  carbonates,  coal,  etc.)^^^ 

'  ™C0  \ 

— I  =  7.23  -  9.66  kg/cm^,  estimated  limits 


'“co. 


8.445  kg/cm^  (computed  from  Ref.  3) 


■“Average 


E.  ELanetaiy  data 


V“e 


r  /r 
®V  ®E 


0.819  (averaged  from  Ref.  3) 

0.950  (estimated  from  Ref.  3,  including  an 
estimate  for  the  height  of  the  cloud  layer) 


®E 


6371  km 


(n) 


sea 
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'Vs 


E 


=  0.90746  (coniputed  from  above) 


980.665  cm/sec' 


,(11) 


F.  Exjjeriment^  (and  highly  uncertain)  data  concerning  the  atmosphere  of 
Venus 


Surface  temperature 

T  =  600°K^^^^ 

s 

Temperattire  of  reflecting  layer  above  -vrfilch  there  are  1CX)0 
* 

m-atm  of  CO^ 

Temperature  at  or  near  the  "cloud  top, ”  but  assumed  to  be 
above  the  300°K  reflecting  layer  (data  compiled  from  original 
references  in  Ref.  6) 


Tg  =  235°K 


THEORETICAL  DEVELOPMENT 
Dole's  Approach 

Dole's  theory  is  based  on  the  assumption  that  since  the  Earth  and  Venus 
are  similar  in  mass  and  size  the  only  differences  betveen  them  must  be  funda¬ 
mentally  due  to  the  amount  of  solar  radiation  received  by  each.  Hence,  an 
evolutionary  theory  of  the  origin  of  the  Earth's  atmosjiiere  expressed  by 
conservation  equations  of  atmospheric  constituents  may  be  looked  at  by  de¬ 
scribing  the  effects  of  placing  Earth  in  Venus'  orbit  some  3  or  4  billion 
years  ago.  Briefly,  hbe  results  are  that  due  to  the  increased  surface  temper¬ 
ature  no  surface  vater  could  exist.  The  water  vapor  would  be  subject  to  a 
_ 

One  m-atm  of  a  gas  is  the  mass  contained  in  a  tube  1  m  long  and  with 
a  1-cmr  cross-sectional  ar«a,  when  the  gas  is  at  standard  pressure  and 
temperature . 
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higher  rate  of  dissociation  due  to  increased  radiation,  and  the  result  would 
be  a  dry  planet .  In  the  absence  of  liquid  water,  CO^  produced  by  volcanic 
action  would  not  go  into  solid  carbonates  or  sea  water.  Atmospheric  turbu¬ 
lence  would  be  very  high,  wind  erosion  would  be  increased,  and,  coupled  with 
a  higher  siirface  temperature  and  normal  isostasy,  surface  oxidation  would  be 
very  great.  Together  with  the  oxidation  of  CO  produced  by  volcanic  action, 
surface  oxidation  would  remove  most  of  the  Og  from  the  atmosphere.  The  bal¬ 
ance  of  Ng  would  be  only  slightly  affected.  Therefore,  the  net  result  is  a 
dry  atmosphere  of  COg  and  Ng.  This  enables  a  quantitative  estimate  of  these 
constituents . 

Under  the  ass-unption  of  static  equilibrtum  with  no  centrtfugal  force, 
and  considering  a  sphertcal  coordinate  system  fixed  to  a  planet,  the  hydro¬ 
static  differential  equation  may  be  shown  to  be  of  complete  generality  as 

dp  =  -  gpdr  (1) 

Now  consider  an  atmosjiiere  of  constant  composition  consisting  of  "i" 
perfect  gases  in  over- all  hydrostatic  equilibrium^  also,  consider  the  ac¬ 
celeration  of  gravity  to  be  constant  through  the  depth  of  the  atmosphere* 

^y  integrating  the  mass  of  the  atmosphere  throxagh  its  depth  above  a  unit 
siirface  area  using  Eq,  (l)  and  perfect  gas  relations,  it  may  be  shown  that 

jmass  of  species  i  in  the  atmosphere!  1  ^i 
\  per  unit  surface  area  1~  g~  14"  ^i 

S  6 

or 


(2) 
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Alternatively,  since 


”l  =  “l^i 


(3) 


Therefore,  If  the  mass  of  species  1  In  the  total  atmosphere  Is  known,  and 
If  the  coD^OBltlon  percentages  are  known  along  with  the  physical  size  of 
the  planet,  the  partial  pressure  of  species  1  at  the  surface  may  be  found. 

Note  that  this  equation  requires  that  each  species  In  a  constant-composition 
atmosphere  be  not  in  hydrostatic  equilibrium  with  Itself,  and  -f^at  the  surface 
partial  pressure  of  a  particular  species  be  not  generally  equal  to  the  wei^t 
per  unit  area  of  the  gas  above  the  surface  as  was  reported  by  Dole.  This  is 
therefore  a  refinement  in  Dole's  findings.  It  is  known  that  for  perfect 
gases  partial-pressure  ratios  represent  ratios  of  percentage  by  volume,  not 
by  weight.  Therefore,  it  would  be  impossible  to  have  the  surface  partial 
pressure  of  a  particular  species  generally  equal  to  the  weight  per  unit  area 
of  this  species  above  the  surface  of  the  planet.  This  refinement  will  de¬ 
crease  the  fraction  of  CO^  reported  by  Dole  by  approximately  4  per  cent. 

The  method  of  computation  then  becomsB 

1.  Consistent  with  the  original  assuii5)ti©n,  the  acceleration  of  gravity 
on  Earth  is  assumed  to  be  a  constant  with  altitude.  From  the  air  composition 

data  on  p.  4  and  from  Eg.  (2),  the  mass  of  N^  in  Earth's  atmosphere  per  unit 
surface  area  is  calculated. 


=  780*11  gm/cm^ 
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2.  From  the  data  on  p.  5  concerning  the  CO^  and  from  step  (l)  above, 
the  coiigxosition  hy  -weight  of  CO^  and  is  found.  Under  Dole's  theory,  gnd 
neglecting  the  trace  of  argon,  this  comprises  -the  atmosphere  of  Venus. 

(%  COg  by  wei^t)  =  100  /  f  +f^] 


(  >  Ng  by  wei^t)  =  100  x  — ^ 


3-  From  step  (2)  abo-ve,  and  the  molectilar-wei^t  data  previously 


mentioned,  the  corr^sition  by  voluiie  is  found. 


(fo  COg  by  -weight)  , 


( CO^  lay  volmae )  =  — - — _ _ 

(%  CO  by  weight) 

1  + - £ - -  - 2_  , 


Ng  by  -volume)  = 


by  weight)  f 


by  wei^t)  f 

^  100  M  ^ 


4.  From  Eq.  (2),  the  molecular  weight  and  planetary  data,  step  (3),  and  the 
assumption  that  the  mass  of  in  the  atmosphere  is  proportional  to  surface  area, 
the  surface  partial  pressure  of  is  found  for  Venus.  Essentially  the  mass 
of  Ng  per  unit  surface  area  is  assumed  to  be  the  seme  for  Earth  and  Venus. 


j  ^  CO  by  -volume 

M  - - _ 

100 


N  by  volume 
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5.  Prom  steps  (3)  and  (4)  the  siirface  partial  pressure  of  cOg  is  found. 

^  COg  hy  volume' 

""COg  "  %  ^  ^  Ng  hy  volume 

'v  ®V 


6.  From  steps  (4)  and  (5)  the  total  stirface  pressure  is  found. 


+  P. 


CO, 


The  resiilts  of  this  computation  are  presented  in  Table  1.  It  should 
not  be  construed,  however,  that  the  significant  figures  shown  are  indicative 
of  the  precision  of  the  theoiy. 


Second  Approach 

For  an  atmosphere  in  adiabatic  equilibrium  the  following  relation  must 

hold 

\  -  s  (*•«) 


■vdiere  a  and  b  denote  any  two  positions  in  the  atmosphere.  Under  our  original 

asstcsrptione  and  with  the  above-mentioned  data,  this  relation  must  hold  below 

the  layer  at  300°K,  above  "vdiich  there  are  ICXX)  m-atm  of  COg.  In  particular, 

it  must  hold  between  the  points  -vdiere  T  =•  300°K  and  T  =  600°K.  Using  the 

J.  8 

data  and  methods  of  Ref.  10,  the  temperatures  T^  and  T  ,  and  Eq.  (4a),  a 

^  s 

relation  among  p^,  p^,  and  the  camposition  is  easily  obtained. 


■vdiere  ^ 


^  dT  +  constant 
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Since  1000  m-atm  of  CO^  =  4.46l  of  CX)^,  tbe  partial  pressure 


of  the  COg  at  point  1  may  be  fotmd  hy  Eq.  (3)  . 


'  n 


CO, 


^CO. 


M  g„ 


V 


Kien  from  the  law  of  partial  preBB^lres,  the  total  pressxire  at  point  1  may 
be  found  as  a  function  of  the  conposition. 


CO, 


+ 


1536  mass  of  in  the  atmosphere  of  Venus  obtained  in  step  (4)  of  Dole's 
approach  is  also  assxaned  to  hold  here: 

I  =  27.845  g-mole/cm^ 

K 


Now,  from  the  law  of  partial  pressures,  Bq.  (3),  and  the  amount  of  N^  the 
total  surface  pressxire  p^  may  be  found  as  a  function  of  composition. 


Hiese  three  relations  yield  a  relation  for  determining  tihe  composition. 


12 


^Jhich  mst  be  solved  for  by  trial.  Canylng  out  the  consultation  gives 
the  following  composition 


COg  a  81.6  per  cent 
Ng  =»  18.4  per  cent 

coBS)ared  to  the  average  values  obtained  by  Dole ' s  approach 

COg  =  87.3  per  cent 

Ng  =  12.7  per  cent 

Hiua  it  is  apparent  that  both  theories  give  quantitative  results  in  close 
a-greement  with  each  other. 

M3DEL  OF  THE  VEMUS  AIMOSFffFRi^ 

Assuming  therefore  that  these  two  theories  do  give  good  quantitative 
estimates  of  the  atmosphere  of  Venus,  the  remaining  problem  is  that  of 
adopting  a  tentative  composition  for  our  model  Venus  atmosphere.  Rather  than 
investigate  the  probable  effects  of  experimental  uncertainty  in  the  data, 
nonunifomn  gravitation,  trace  elements,  impeirfect 


gases,  etc.,  idiich  would 
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in  fact  invalidate  all  of  the  original  assurt^tions  and  require  more  accuracy 
than  is  obviously  contained  in  this  model,  ve  make  a  fundamental  assumption: 
that  the  temperature  lapse  rate  with  altitude  is  slightly  less  than  dry  adi¬ 
abatic.  The  experimental  data  and  all  of  the  original  assumptions  are  re¬ 
tained  except  that  of  a  dry  adiabatic  lapse  rate  below  T^.  As  an  alterna¬ 
tive  to  Eq.  (l+a),  the  well-known  relation 


could  have  been  written  for  an  adiabatic  atmosphere  of  perfect  gases.  It 
^vas  not  used  because  y  has  a  strong  variation  with  temperature  for  this 
mixture,  and  a  "proper"  average  y  would  have  had  to  be  chosen  as  in  the  "C" 
data  on  p,  5.  if  a  polytropic  atmosphere  is  assumed  where  an  adiabatic  atmos. 
phere  was  previously  assumed,  the  following  equation  holds  by  definition 


Pb 


(5) 


For  an  adiabatic  atmosphere,  q  =  7  ;  for  an  isothermal  atmosphere,  q  =  1,- 
therefore,  we  will  choose  1  <  q  <  7.  Using  this  criterion  and  Eq.  (5)  m 
place  of  Eq.  (1+a),  and  repeating  the  computations  of  the  second  approach, 
it  is  found  that  the  conqjositlon  must  contain  more  than  8I.6  per  cent  CO 
by  volume.  This  tends  to  reduce  the  discrepancy  between  the  two  approaches. 
Ihe  atmosplriere  of  Venus  is  now  arbitrarily  chosen  to  consist  of  85  per  cent 
COg  and  I5  per  cent  by  vol-une. 

With  a  composition  assumed,  the  atmosphere  is  modeled  as  follows: 


Il4- 


1.  Fram  the  known  compoeltlon,  aiaount  ©f  and  Eq..  (3),  the 
total  surface  pressure  is  found. 

2.  From  the  known  eompositian,  amount  of  GO^  above  point 

and  Eq.s.  (3)  and  (5),  the  exponent  q  is  fo\md. 

3.  This  polytropic  atmosphere  is  assinued  to  hold  true  until  point 
2  is  reached,  -vdiere  T^  =  235°K,  the  ass\imed  "cloud  top." 

h.  The  atmosphere  will  be  considered  isothermal  above  point  2, 

■vdiere  1=1^  =  235°K  and  q  =  1.  Now  it  is  highly  improbable 
that  this  situation  precisely  prevails,  but  it  should  represent  a 
reasonable  assun^rtion  over  altitudes  of  aerodynamic  interest.  From 
the  standpoint  of  radiative  equilibrium,  the  temperature  probably 
drops  to  about  190°K  at*  55  km.  There  may  also  be  a  sli^t  increase 
in  temperature  above  this  level  due  to  0^  heating.  This  will 
be  much  more  suppressed,  however,  than  the  case  on  Elarth  because 
of  the  absence  of  much  0^  at  these  altitudes  (provided,  of  course, 
that  there  is  no  appreciable  amount  of  Og  in  the  atmosphere  of 
Venus) .  Effects  of  this  kind  are  neglected  in  this  model,  how¬ 
ever,  because  they  are  inconsistent  with  the  accuracy  possible 
through  use  of  the  other  assumptions  of  the  model. 

5.  With  the  above  inf9rmati*n,  Eq.  (5),  and  the  perfect  gas  law, 

Eq.  (1)  can  be  integrated.  This  yields  a  complete  pressure-densily- 
ten^rature-sQ-titude  relation  consistent  with  the  assumptions  and 
data.  This  relation  is  presented  below,  and  it  is  represented  in 
Table  2  and  in  Fig.  1. 


Venus 
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Tdiere  consistent  with  the  assTiJi5)tionB  and  derivation  of  the  model 

T  =  600°K 
s 

=s  6,78^  atm 

Pn  =  5-73^  X  10"^  g/cm^ 

r  =  6052.45  km 
s 

Pg  =  5*^*^  X  10”^  atm 
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Pg  =  1.179  X  10"^  g/cm^ 

gg  =  889.9  cm/sec^ 

=  5.1443 

M  =  41.612  g/g-mole 

=  r  +  42.16  km 

2  s 

R  =  8.31433  X  10^  erg/g-mole  °K 
DISCUSSION 

It  is  necsessary  at  this  point  to  review  exactly  lAiat  has  been  done  in 
order  to  see  idiat  assuz^ions  are  basic  to  this  atmospheric  model.  In  analogy 
to  the  Earth,  it  is  immediately  evident  that  the  assun^rtions  of  nnifonn 
gravitation,  constant  composition,  perfect  gases,  and  hydrostatic  equilibrium 
with  no  centrifugal  force  are  not  the  major  sources  of  uncertainty,  it  is 
in  the  statements  that  followed  that  great  uncertainty  was  introduced.  In 
essence.  Dole's  theory  was  assumed  correct  at  the  outset,  and  was  investi¬ 
gated  as  to  compatibility  with  some  experimental  observations.  It  was  neces¬ 
sary  at  the  start  to  assume  that  CO^  and  N^  are  the  only  major  constituents 
in  the  atmosphere  and  additionally  to  assure  the  total  armrunt  (not  proportion) 
of  in  the  atmosphere.  Then  it  was  necessary  to  assume  accuracy  in  sque 
er^rimental  data  that  are  known  to  contain  a  hi^  degree  of  uncertainly. 
Further,  it  was  necessary  to  speculate  on  the  nature  of  the  cloud  layer  sur¬ 
rounding  the  planet  so  that  the  nature  of  the  temperature  lapse  rate  with 
altitude  could  be  specified.  At  this  point  it  would  appear  that  the  entire 
problem  had  been  assumed  away.  However,  there  are  indications  that  this 
not  be  the  case. 


may 
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Consider  first  that  Bole’s  approach  is  accepted  and  that  the  question 
is  one  of  corapatihili-t^  with  some  experimental  measurements.  In  dealing 
with  the  second  approach,  it  is  found  first  of  « l i  that  any  uncertainty  in 
the  amount  of  chosen  greatly  affects  the  derived  composition.  Uncertain¬ 
ties  of  1  per  cent  in  this  amount  affect  the  resulting  fraction  of  by 
approximately  1-2  per  cent.  It  is  also  found  that  uncertainties  in  the 
data  concerning  the  amount  of  CO^  above  the  "reflecting  surface"  greatly 
affect  the  derived  composition.  For  Instance,  if  there  were  only  5OO  m-atm 
©f  GQg  above  this  svirfaee  Instead  of  the  chosen  flgme  of  1000  m-atm,  the 
derived  composition  wotild  contain  approximately  3O  pesr  cent  by  volume. 

Die  results  are  also  highly  dependent  upon  the  ter^ratures  chosen  and  i^pon 
the  nature  of  tei^rature  lai>se  rate  below  the  "reflecting  surface.”  In 
other  words,  there  are  no  chosen  parameters  that  tend  to  suppress  uncertainly 
effects  or  that  cause  less  uncertainty  in  the  results  for  composition 
than  is  contained  in  the  chosen  parameter.  It  is  then  apparent  that  the 
close  agreement  of  the  two  approaches  must  either  be  a  remarkable  coincidence 
or  have  some  basis  in  fact.  However,  this  model  is  inconsistent  vlth  the 
measurements  of  Ref.  "J,  and  consequently  with  the  results  of  Ref.  8. 

It  must  also  be  noted  that  the  amGunt  of  reported  by  Dole  is  approxi¬ 
mately  one-half  that  which  would  result  if  the  figures  for  the  amount  of 
Earth  00^  from  Ref.  14  were  used.  Attention  must  also  be  called  to  Ref.  I5 
in  Tdileh  a  method  similar  to  the  second  approach  here  is  used  to  derive  a 
coinpositlon  of  approximately  75  per  cent  COg  and  25  per  cent  Ng. 

All  that  is  claimed  here  is  that  under  the  assuiuptlens  a  tentative 
atmosphere  of  Venus  has  been  derived.  Now  if  the  asstm^ptlons  involved  tiltl- 
mately  prove  reasonable  and  the  ejqperimental  data,  and  consequently  the 
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conclusions  derived  therefrom,  are  reasonably  correct,  then  the  close  agree- 
neat  of  the  tw  asproa^s  Inileates  that  the  <*oeen  caopoeltlon  ^  he 
representative  of  the  Venus  atoosphere.  Further,  the  pressure>tsnpisrature 
distrihutlon  helov  the  cloud  layer  may  be  reasonably  oorreet.  Oie  diffi¬ 
culties  in  accepting  an  isothermal  atiaosphe:re  above  the  cloud  layer  have 
been  diseussed  previously.  Above  100  tan,  photodissociation  and  ionisation 
will  surely  beeeme  appreciable,  and  in  going  to  these  higher  altitudes  the 
assTssptlon  of  a  uniform  gravitational  field  becomes  poorer.  Ube  cutoff  at 
100  tan  was  chosen  because  the  density  at  this  point  roughly  corresponds 
to  that  density  in  the  Earth's  atmosphere  (at  approximately  300,000  ft) 

at  which  incipient  aerodynamic  effects  associated  with  re-entry  nose  cones 
are  no  longer  negligible. 

One  immediate  objection  to  the  acceptance  of  this  model  concerns  the 
arbitrary  choice  of  the  teanperature  lapse  rate  below  the  cloud  layer.  In 
reality,  the  real  Justification  for  choosing  a  lapse  rate  slightly  less  than 
an  adiabatic  lapse  rate  Is  that  this  choice  reduces  the  discrepancy  in  atmos¬ 
pheric  eati5>o3ltlon  between  the  two  approashes.  However,  it  may  be  partially 
Justified  by  other  arguments.  For  instance,  if  the  cJ.®uds  are  composed  of 
water  and  ice  crystals,  the  lapse  rate  in  the  cloud  layer  must  be  less  than 
the  dry  adiabatic  lapse  rate  becaorse  of  the  release  of  the  heats  of  vapori¬ 
zation  and  fusion  as  altitude  Increases.  The  chosen  lapse  rate  may  then  be 
considered  an  average  between  the  sxrrface  and  the  top  of  the  cloud  layer. 

Yet  another  possibility  may  be  considered.  If  the  surface  temperature  is 
subtly  greater  than  600°K,  as  may  well  be  possible  considering  the  un¬ 
certainty  in  the  data.  Eg.,  (l»a)  requires  that  the  surface  pressure  be  higher 
than  that  originally  derived.  Assuming  that  the  amount  of  is  fixed,  Eq.  (a) 
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xequlzna  that  there  he  mere  CX)^.  !Ihe  ohosea  lepae  rate  nagr  then  he  een- 
sldered  as  a  eexreetlen  device  to  fix  the  surfaee  teagerature  at  6O0®K,  but 
still  he  censlsteat  -with  -fche  amount  of  CO^  In  the  atnosjhere.  Of  course 
this  argument  vezicB  in  reverse  If  the  surface  -benyerature  Is  less  -than  600®K, 

It  dx>e8  not  seem  profltehle  -to  investigate  all  of  these  minor  effects  in  -view 
of  the  uncertainties  already  introduced.  Oils  lapse  rate  -was  chosen  singly 
because  it  serves  -to  reduce  -the  discrepancy  between  -the  two  approaches  aoid 
appears  -to  he  as  good  as  any  other  assun^ptlon  at  -this  tine. 

It  is  Interesting  -to  note  -that  If  the  clouds  surrounding  Venus  are  as¬ 
sumed  "to  consist  of  lee  crys-tsils,  according  to  "the  present  model  -the 

vapor  p3^ssure  above  lee  at  »«  235°K  indlea-tes  -that  -vrater  -yaper  comprises 
0.29  per  cent  by  -volume  of  -the  Veniis  a-tmosihere  abo-ve  -the  cloud  layer  (as¬ 
suming  a  cons-tant  mixing  ratio  abo-ve  -this  point).  This  la  larger  by  a  factor 
of  approximately  tO  -than  that  -vdilch  -would  exist  in  -this  model  according  to 
the  recent  measurements  of  -the  Moore-Ross  balloon  flight. It  must  be  cau¬ 
tioned,  however,  that  -these  meastirements  contained  a  high  degree  of  uncertain-ty. 

In  view  of  -the  above-mentioned  difficulties,  -this  model  is  proposed  only 
as  a  first  approximation  -wl-th  vdilch  aerodynamic  entry  problems  concerning 
Venus  may  he  studied.  It  may  be  used  as  a  mean  to  which  pert-arbatlons 
may  be  applied;  this  mean  may  be  bracke-ted  by  ©-ther  existing  theories  until 
more  experimental  data  are  accumulated.  Ihe  significant  figures  presented 
for  this  model  are  ob-vlously  not  indlcatl-ve  of  the  precision  of  -the  Isnowledge 
©f  -the  Venus  atmoBihere,  but  are  consis-tent  -wl-th  -the  proposed  model  if  -the 
experimental  data  and  assumptions  in-vol-ved  are  -taken  as  "exact." 
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III,  TEIE  IHERMDDYNAMIC  PROPERTIES  QF  15  PER  CENT  N 

Airo  85  PER  CENT  CO„  BY  VOLUME 

"""  "  ' 

CONDITIONS  AMD  ASSUMPTIONS 

The  thermodynamic  properties  of  a  mlxtiare  of  85  per  cent  CO^  and  I5  per 
cent  Ng  by  volume  will  be  presented  on  the  premise  that  they  may  be  useful 
regardless  of  specific  application  to  this  tentative  Venus  atmosphere.  In 
this  light,  the  calculations  of  this  section  will  be  carried  out  for  exactly 
the  stated  composition  and  will  contain  far  less  computational  imcertalnty  than 
does  Section  II  concerning  the  constituency  of  the  atmosphere  of  Venus. 

The  gas  mixture  constituents  are  assumed  to  obey  the  perfect  gas  lav, 
pM/RpT  —  1.  Investigation  of  the  critical  constants  of  the  constituents 
shows  that  this  is  a  good  assumption  over  the  selected  teirperature  range  of 
from  1000  K  to  2^,000*^.  ^^2  most  from  the  perfect  gas  law  (the  ef¬ 

fects  of  compressibility  for  this  constituent  have  been  discussed  by  Ray- 
mond) .  Since  is  better-behaved  in  this  respect  than  CO^,  the  com¬ 

pressibility  at  1000°K  and  10^  atm  for  this  mixture  is  slightly  depressed 
from  the  value  of  1.026  quoted  by  Raymond  for  pure  COg,  However,  this  as- 
sumptlon  becomes  poorer  at  temperatures  below  1000°K  and  at  the  higher  pres¬ 
sures.  The  Justification  for  neglecting  the  effects  of  intermole cular  forces 
at  these  lower  temperatures  is  that  the  thermodynamic  states  of  interest  with 
respect  to  the  atmosphere  of  Venus  do  not  occur  where  large  effects  of  Inter- 
molec\alar  forces  are  apparent.  When  the  pressure  is  nearest  the  critical 
pressure  (for  COg)  the  ten5)erature  is  high;  where  the  temperature  approaches 
and  falls  below  the  critical  temperature,  the  pressure  is  low.  Again,  the 
presence  of  Ng  will  generally  tend  to  influence  the  mixture  to  behave 
like  a  perfect  gas  than  will  pure  COg, 
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The  mixtiire  will  be  assuned  -fco  be  In  thermodynamic  equilibrium.  From 
the  standpoint  of  the  use  of  these  properties,  this  assumes  that  for  aigr 
particular  state  to  \dilch  these  properties  are  to  be  applied,  e.g.,  aft  of 
a  shock  wave,  the  time  reqioired  for  relaxation  phenomena  to  take  place  is 
much  shoiter  than  that  in  \diich  the  gas  changes  appreciably  in  state.  The 
equilibrium  catjvsltion  of  the  gas  mixture  is  then  computed  by  the  method 
of  Ref.  18  under  the  criteria  that  the  Gibbs  Free  Energy  is  a  minimum  for  a 
selected  particular  temperature  and  pressure. 

Below  1000°K  the  mixture  is  assumed  to  consist  solely  of  molecular  CO 

2 

and  N^.  This  ass\imption  is  Justified  by  noting  the  calculated  equilibrium 
con^-T'fiitlon  at  T  =  1000*^K  and  p  =  10  ^  atm.  Using  perfect  gas  relations, 
the  thermodynamic  properties  of  the  mixture  can  be  calculated  from  the  basic 
data  in  Pief.  10.  However,  a  slight  correction  must  be  applied  to  these  data 
in  order  that  they  may  become  consistent  with  the  data  above  1000°K;  the 
method  of  correction  will  be  described  later. 

METHOD  OF  CALCULATION 
T  >  1000°K 

Under  the  criteria  of  significantly  contributing  to  the  thermodynamic 
properties  of  the  mixture  it  was  first  as sinned  that  the  mixture  could  con¬ 
sist  of  the  following  possible  species:  00^,  0^,  NO^,  CO,  NO,  0^,  C,  N, 

0,  C  ,  N  ,  0  ,  Og,  0  ,  CO  ,  NO  ,  Ng,  Og,  C  ,  N  ,  O"*^^,  and  e  ,  Using  the 
basic  referenced  data  in  Tables  3-6,  and  the  n^thod  of  Ref.  18,  the  equi¬ 
librium  congiosition  was  calculated  over  the  teii5>erature  range  1000°K  <  T  < 
2i(-,000°K.  Constituents  present  in  amounts  smaller  than  10~*^  g-mole  per  orlgi' 
nal  "cold"  g-mole  of  mixture  were  then  considered  negligible  end  not  listed. 
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Thus,  thermodyBaiiic  properties  are  presented  for  the  foUevlni;  significant 
censtltuents  In  the  temperature  ranges  specified  below: 

1000®K  <  T  <  li.OOO®K 


V  V  ° 

4000°K  <  T  <  24,000°K 

CO^,  0^,  CO,  NO,  C,  N,  0,  n'*',  0^ 

0“,  CO^  HO"^,  0*,  C"^,  n"^,  0“^,  and  e" 


Uie  eqTollibrlum  conrposition  is  presented  on  the  basis  of  one  g-mole  of  the 
cold  mixture  in  Figs.  2  and  3  snd  Tables  7  and  8  over  a  temperature  range 
of  1000°K  to  24,000°K  and  a  pressure  range  of  10^  atm  to  10~^  atm. 

Since  the  thenaodynamic  data  are  presented  on  the  basis  of  one  g-mole 
of  the  "cold"  gas  mixture,  i.e.,  14-1.612  g,  the  following  relation  applies 


since  n  -  1  and  £  n  -  n.  Thus,  n  /n  represents  the  mole  fraction  of 
i  ^ 

species  i.  The  datum  density,  p^,  is  chosen  to  correspond  to  p^  =  1  atm  (on 
the  surface  of  the  earth)  and  T^  =  273.l6°K  for  the  original  "cold"  mixture. 
Thus 

p  M 

Pq  ~  ~  1.85647  X  10  ^  gm/cm^ 

o 


TOien 


£_ 

p 

o 


(p/Po)(M/Mq) 


(12) 


24 


Fig.  2  —  Total  number  of  moles  as  a  function  of  temperature 


Tempera 


27 


Temperature  ,  T  (thousands  of  °K  ) 
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10 


20 


Temperature,  T  (thousands  of  °K) 


Fig.  3e  —  Molar  composition 
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Fig.3f — Molar  composition 


Fig.  3g —  Molar  composition 
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The  molecular  vei^ts  and  densities  are  shown  in  Tables  9  and  10  and  Figs. 

4  and  5« 

The  datim  state  chosen  for  the  internal  energy  equal  to  0  is  where 
carbon  (graphite),  gaseo\rs  molecular  oxygen,  and  gaseoiis  molec\ilar  nitrogen 
exist  at  a  temperature  of  0°K  and  zero  pressure.  Then  the  internal  energy 
per  "cold"  g-mole  of  the  mlxtiire  may  be  calculated  as  follows: 

E  =  £  n  E  (13) 

i  ^ 

vbere 

E^  =  (E°  -  E°)^  +  E°^  (14) 

and  E°  is  the  heat  of  formation  of  species  i  from  the  grotind  state  mentioned 
above.  In  dimensionless  form,  the  computing  equation  becomes 


The  basic  data  are  taken  from  the  previously  mentioned  tables. 

Since  the  enthalpy  is  given  by 

Hi  =  E^  +  p^V  =  E^  +  n^  RT  (l6) 

the  computing  equation  is 


20 


25 


10  15 

Temperature,  T  (thousands  of  °K) 

—  Molecular  weight  versus  temperature 
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The  Gibhs  Free  Energy  is  given  by 


Also 


Then 


£n  +  1 


or 


(18) 

(19) 

(20) 


Since 


and 


the  computing  equation  is 


(21) 


(22) 


(23) 


Again,  the  basic  energy  data  are  taken  from  the  previously  mentioned  tables. 

The  values  for  internal  energy,  enthalpy,  and  entropy  are  sho’.'m  in 
Tables  11-13  and-  Figs.  6-8,  respectively.  Hi^-temperature  Mollier 
charts  are  presented  in  Figs.  9  and  10. 


6500 


<  18.000 


Enthalpy,  H/RTo 


^UUU 


Entropy  ,  S/R 


Fig.  10  Mollier  chart  to  12,000  °K  for  a  volumetric 
mixture  of  85%  CO,  and  15%  N, 


kk 


To  convert  tlie  ttienoodynaicic  propei^iies  to  the  mass  basis^  it  is  inerely 
necessary  to  divide  the  properties  by  the  "cold"  molecular  weight.  Thus 


_e _ 1^ 

RT  "  RT„ 
o  o  o 


(25) 


h  1_  H 

HP  “  M  RT 
o  o  o 


and 


s  S 

RT  “  M  R 
o  o 


(26) 


(27) 


T  <  1000°K 

The  thermodynamic  properties  in  this  range  were  con^juted  from  Tables 
12  and  17  of  Ref.  10  for  a  mixture  of  O.85  g-mole  of  CO^  and  0.15  g-mole  of 
Ng.  In  this  reference  the  zero  for  internal  energy  and  enthalpy  is  taken 
at  0°K.  It  was  necessary,  therefore,  to  add  the  heat  of  formation  of  COg 
to  the  enthalpy  -of  COg  tabulated  in  Ref.  10.  In  order  for  these  data  also 
to  match  the  coniputed  properties  at  1000°K  but  still  retain  the  zero  at  0°K, 
it  is  necessary  to  multiply  the  enthalpy  data  in  Ref.  10  by  a  factor  only  very 
slightly  different  from  unity.  In  effect,  this  means  that  the  specific  heat 
at  constant  pressure  is  multiplied  by  this  factor.  The  corapiuting  equation 
then  becomes 


RT  RT 
o  o 


^CO  ^CO 

0.85  1  — -  -  173-1048 


+  0.15 


K^h. 


2 _ 2 

1.8  , 


(28) 


where  the  K's  are  the  multiplying  factors  to  match  the  data  at  1000  K  and 


=  1.000324 

le  =  1.0004675 
2 


^5 


Now  the  entroEjr  function  ®  ,  tabulated  In  Bef.  10,  represents  the 


following: 


E°  -  F°  +  R 


.  -  f 


c 

^  dT  +  constant 


Thus,  the  same  multiplying  factor  must  be  applied  to  this  function  as  Is 
applied  to  the  enthalpy  above.  The  final  matching  of  these  data  with  the 
computed  data  at  1000°K  Is  then  taken  care  of  by  the  additive  constant  in 
(30)-  Ihus,  the  con^niting  equation  is 


1  =  0. 


^co  ®co. 


''''2  2  «  N 

•  ^  R  *  ®C0  ■  (0-85p)  +  0.15  — ^  +  Bjj  -  in  (O.lJ+p) 


■idle  re 


=  0.0114 


Sfj  =  -0.0008 
2 


The  values  for  the  enthaljy,  internal  energy,  and  entropy  over  the  range 
of  150  K  to  800°K  and  10^  to  lo"^  atm  are  shown  In  Figs.  6  -  8  and  Tables  11  -  13. 
In  Table  2,  values  are  also  included  for  the  entropy  at  the  particular  pressure- 
temperature  combinations  that  occur  in  the  model  of  the  Venus  atmosphere  se¬ 
lected  in  Section  II. 


DISCUSSION 

Since  much  of  the  basic  data  were  taken  from  Gilmore’s  paper, the 
reader  is  referred  to  Ref.  9  for  an  adequate  discussion  of  the  accuracy  of 
the  basic  data.  Constituents  having  a  large  error  in  the  basic  data  (due 


k6 


•to  a  siiramatlon  o-ver  too  fev  energy  le'vels  in  the  in-ternal  partition  ftmction) 
are  present  in  such  small  amounts  at  "the  ■temperatures  and  pressxjres  -^daen  -the 
error  is  quite  significant  -that  -the  error  introduced  in  -the  -thermodynamic 
functions  may  be  considei-ed  -veiy  slight.  It  shoxold  also  be  mentioned  -that 
■the  data  are  significantly  in  error  near  -the  liquefaction  points  of  -the  COg. 

For  more  accura-te  results  in  this  region,  bet-ter  data  for  CO^  should  be 


used. 


Selec'ted  cases  rere  calcula-ted  oi/er  -the  coraple-te  pressure  range  at  -temper¬ 
atures  of  5000°K  and  12, 000°K  for  a  mixture  of  84.7  per  cent  CO^,  I5  per  cent 
Ng,  and  0. 3  per  cent  HgO  vapor  by  volume.  This  was  done  to  investigate  the 
effect  of  wa-ter  -vapor  on  -these  properties  in  -the  event  1:hat  a  substantial 
amount  of  -wa-ter  -vapor  actiiaJly  exists  in  "the  a-tmosphere  of  Venus.  The  ad¬ 
ditional  constituents  considered  -were  H^O,  H^,  OH,  H,  H^,  OH''',  H",  and 
OH  .  Table  8H  sho-ws  the  effect  of  -this  addition  on  -the  composition  and  -thermo¬ 
dynamic  properties  of  the  mixtiire.  It  may  be  s-ta-ted  that  -the  -thermodynamic 
properties  are  not  affected  by  more  -than  -the  order  of  percen-tage  addition  of 
H^O. 

The  three  most  important  differences  between  this  case  and  that  of  pure 
( 17 1 

COg^  '  are  the  reduction  in  internal  degrees  of  freedom  at  low  temperatures, 
since  Ng  is  present;  the  appearance  of  NO  at  moderate  temperatures;  and  the 
reduction  in  electron  concentration  at  high  temperatures  due  to  the  high  ion¬ 
ization  potential  of  Ng.  Since  NO  forms  quite  readily  around  5000°K,  more 
energy  goes  to-ward  dissociation  at  this  level  than  in  the  pure-COg  case. 

Since  NO  ionizes  quite  readily,  the  electron  concentration  at  this  temperature 
level  is  higher  than  that  for  pure  COg .  All  effects  mentioned  are  generally 
of  the  order  of  the  percentage  of  Ng  addition  to  pure  COg. 
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IZi _ NORMAI^  SHOCK- WAVE  OLARACTBRIgTICS  OF  A  TBUTATIVE  VEHUS  ATMOSPHERE 

coiroinoHS  and  assumptions 

The  gasdynamic  properties  of  the  atmosphere  of  Venus  depend  upon  the 
gas  mlxt\ire  assumed  and  the  thermodynamic  properties  at  the  hi^  temperature 
and  pressure  ratios  usually  ejqpected  throx;igh  a  normal  shock  at  high  entry- 
speeds.  Therefore,  it  should  be  recalled  that  the  gas  mixture  assrnned  is 
85  per  cent  CO^  and  I5  per  cent  by  volume  (Section  II),  and  that  -the 
calculations  of  the  thermodynamic  properties  are  based  upon  the  assumption 
of  -tbermodynamlc  equilibrium  (Section  III), 

In  addition,  it  must  be  assumed  that  aerodynamics  of  continuum  flow 
applyj  that  is,  -that  -the  shock  thickness  must  be  very  much  smaller  thain  the 
charac-teristic  length  of  the  object  causing  -the  shock.  Ihis  is  necessary 
in  order  -that  -the  conservation  equations  may  be  applied  across  a  -very  thin 
flow  discontinuity.  Since  -the  highest  altitude  (and  lo-west  density)  chosen 
corresponds  -to  -the  altitude  in  Earth’s  atmosphere  -vdiere  this  assumption  still 
generally  applies,  the  assumption  does  not  impose  an  appreciable  limitation. 

However,  it  must  be  cautioned  that  the  mean  free  path  at  100  km  is  of  the 
order  of  10  cm . 

METHOD  OF  COMPUTATION 

Several  methods  of  calculating  shock-wave  characteristics  were  investi¬ 
gated  •  ^  9,  20,  21  )  ,  _ 

ganea,  however,  because  the  upstream  state  must  be  specified  and  be¬ 

cause  of  the  desirability  of  a  rapid  convergence  of  the  iteration  procedure, 
me-thod  of  Ref.  2O  was  selected,  since  it  api>eared  -to  be  well-suited  to  the 
solution  of  the  problem.  Applying  conservation  of  mass,  momentum,  and  energy 
across  a  flow  discontinul-ty 


Continul-ty  ; 


^1^  =  ^2^2 


(32) 


Mooeatum: 


-  Pg  +  PgUg 


(33) 


Energy: 


1  ..2 


-"a"! 


B  +  - 


1  ..2 


*2  2^ 

Siibstituting  Eg.  (32)  in  Eg.  (33)  and  jnanipulatlng,  the  presaure  ratio 
becomes 


(34) 


^2 


“o 


,  »->  2  I  ,  ^ 


Rearranging  and  manipulatiiig  Eg.  (34),  the  enthaliy  becoines 


H, 


KT  "  KT  RT  2 
000 


1  - 


(35) 


(36) 


Fram  Section  III  the  eguation  of  state  is  presented  in  Tables  10  and  12  in 
the  form 


RT  -  ^  (^2^  ^2/^0^ 


(3Ta) 


P 

P 


2 

o 


f  (V  Pg/p^) 


(37b) 


Finally,  from  Eg.  (32) 


J^9 

With  the  equatlms  above,  the  machine  Iteration  procedure  foUovm  this 
pattern: 

1.  FlcJc  an  upatream  thermodynamic  state  and  veloci-^  u^,  and  entei* 
all  known  quantities... 

2.  Aasuns  '‘12/'^  “  Calculate  Pg/p^  from  Eq.  (35)  and  Hg/RT^  from 
Sq.  (36). 

3.  Enter  these  values  in  the  tabular  relations  (Eq.  (37))  and  find 
^2/^0- 

k.  Calcxilate  u^/u^  frem  Eq.  (38)  and  enter  this  into  Eqs.  (35)  sind  (36). 

5.  Repeat  this  procedure  imtil  u^/u^  obtained  in  step  (k)  above  differs 
from  that  \jsed  in  the  prevlatis  iteration  by  less  than  O.OCOl. 

6.  Print  out  Pg/Pj.^  '^2^'^r  ^2!^' 

The  results  of  these  computations  are  presented  as  a  function  of  altitude 
in  the  adopted  model  of  the  Venus  atmosphere  in  Figs.  11a  -  lid  Tables 
14a  -  l4d  over  a  flight-speed  range  of  2,000  to  40,000  ft/sec.  A  typical 
sound  speed  is  9IC  ft/sec  (at  an  altitude  of  34.65  km  in  this  model)  so 
that  the  free- stream  Ifech  number  range  is  ro\:ghly  2.2  to  44. 

DISCUSSION 

It  is  interesting  to  note  that  to  a  good  approximation  at  low  Mach 
numbers,  air  may  be  treated  as  a  perfect  gas  with  constant  specific  heats 
as  far  as  shock-wave  relations  are  concerned;  however,  this  atmosphere 
cannot  be  so  treated  because  of  the  marked  variation  with  temperatixre  of 
the  specific  heats  of  CO^  in  the  normal  temperature  range  of  interest.  It 
may  be  shown  throu^  an  error  analysis  that  the  iteration  procedure  can  be 
cut  off  -vdien  two  successive  values  of  u^/u^  differ  by  less  than  0.0001  muS 
yet  generally  yield  errors  in  the  thermodynamic- state  variables  of  the  oarder 
of  0.01  per  cent. 


Pressure  ratio,  p^/ 
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Fig.l  I  a -2— Normal -shock- wave  characteristics 
of  a  tentative  Venus  atmosphere 
(pressure  ratio,  P2/P1) 


Temperature  ratio, 
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Fig.  Mb  —  Normal-shock-wave  chorocteristics  ofo  tentative 
Venus  atmosphere 
(temperoture  ratio,  Tg  /T|) 


fig.  lid  I  Normol-shock-wove  chorocteristics  of  o 
tentative  Venus  atmosphere 
(velocity  ratio,  ug/u,  ;  U|,from  0  to  40,000  fps) 


APPENDIX  A 
TABLES 
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Average  8.1022  6157.7  100.000  I  100.000 
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Table  3 

BASIC  PHYSICAL  CONSTANTS  AND  CONVERSION  FACTORS 

Universal  Gas  Constant  (chemical  scale) 

R  =  8.31433  X  10^  erg/g-imole  °K 
=»  1.98717  cal/g-mole  °K 
«  82.0561  atm  cm^/g-mole  °K 

Avogadro's  Number  (chemical  scale) 

L  =  6.02306  X  10^~'/g-mole 


Energy  Conversion  Factors  (chemical  scale) 

1  cal/g-mole  •»  4.33605  x  lO”^  ev 
1  thermochemical  calorie  ■  4.18400  x  10^  erg 

Pressure  Conversion  Factor 

6  P 

1  standard  atmosphere  ■  1.01 325  x  10  dyne /cm 

Ratio  of  Physical  to  Ghemioal  Scales 
^jhys^chem  "*  ^phys^chem 


1.000275 
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Table  4 


ENEilGIES  OF  FORMATION  E° 

o 


Constituent 

0 

kcal/g-mole 

All  obtained  from  Ref.  9 
by  method  indicated  below 

“2 

-93.9639 

Direct 

^2 

0 

Definition 

O2 

0 

Definition 

CO 

-27.1992 

Direct 

NO 

21.479 

Direct 

C(g) 

169.99 

Direct 

N 

112.535 

Computed 

0 

58.985 

Computed 

C^ 

429.832 

Computed 

N-" 

448.051 

Comimited 

-  + 

0 

373.033 

Computed 

_  + 

°2 

277.9 

Computed 

co'*' 

295.981 

Computed 

NO 

234.879 

Computed 

++ 

C 

992.125 

Computed 

N 

1130.969 

Computed 

„++ 

0 

1183.77 

Computed 

0" 

25.485 

Coniputed 

e" 

0 

Definition 

denotei  extrapolation  or  interpolation  Veo  from  r.  R.  Giljnore,  unpublished  thensodynamic  properties  from  statlitlcsj.  mechanic: 
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Table  8c 

lOLAH  OOMPOSITIOS  n,  (PRESSURE  -  1.0  AM) 
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70 


72 
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Table  8l 


THERMODYNAMIC  PROPERTIES  FOR  ATMOSPHERE 
WITH  0.3  PER  CENT  WATER  VAPOR 


Table  8l 


thermodynamic  properties  for  atmosphere 

WITH  0.3  PER  CENT  WATER  VAPOR 


Property 


M 

P/Po 

e/rt 

0 

T  =  5000°K 

16.211 

19.494 

21.438 

22.135 

22.529 

23.594 

26.752 

2.1323"^ 

2. 5642 

2.8193"^ 

2.9115"^ 

2.9634“^ 

3.1035"^ 

3.5189 

459.25 

257.77 

170.54 

145.09 

134.78 
117.16 
74.32 

506.15 

296.76 

206.00 

179.44 

168.52 

149.39 

102.74 

s/r 


96.56; 

79.7^1 

70.14i 

6^.317 

59.^37 

5^.222 

47.8J12 


T  =  12000°K 


7.2981 

7.^171 

8.2261 

10.319 

12.10^2 

13.784 

14.642 


3.9999"^ 

2469.4 

2739.4 

172.87 

4.0651"° 

2429.4 

2675.4 

158.38 

4.5084"5 

2072.7 

2294.5 

137.31 

5.6554"^ 

1420.0 

1596.8 

110.99 

6.8l88"3 

1006.8 

U-53.5 

92.503 

7. 5544 

822.47 

954.85 

80.722 

8.0246"^ 

731.09 

855.70 

71.731 

-«  CO  in  CVJ  w  -I  -t  rA  kA  cO  cO  OJ* 


Belov  lOOO'K  this  relation  Is  given  by  p/p  =  275.16 
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Table  Il)-A 

NORMAL  SHOCK  WAVE  CHARACTERISTICS  OF  A  TENTATIVE  VENUS  ATMOSPHERE 

PRESSURE  RATIO,  Pg/Pj 


Velocity 

1  Altitude  ( Kin^ 

(ft/sec) 

0 

23.10 

42.16 

60 

80 

100 

20  X  10^ 

2.692 

4.016 

6.773 

6.775 

6.773 

25  X  lo; 

4.275 

6.374 

10.745 

10.749 

10.745 

30  X  lOp 

6.217 

9.269 

15.638 

15.644 

15.637 

35  X  lOp 

8.504 

12.709 

21.457 

21.466 

21.456 

1*0  X  lOp 

11.132 

16.669 

28.208 

28.219 

28.206 

1*5  X  lOp 

14.135 

21.136 

35.838 

35.853 

35.840 

50  X  lOp 

26.178 

44.343 

44.365 

44.359 

55  X  lo; 

31.840 

53.921 

53.951 

53.952 

6o  X  lOp 

37.932 

64.381 

64.471 

64.526 

65  X  lOp 

44.687 

75.769 

75.943 

76.074 

70  X  lOp 

51.941 

88.288 

88.607 

88.799 

75  X  lo; 

59.702 

101. 501 

101.893 

102.744 

80  X  lOp 

68.066 

115.679 

115.945 

117.192 

116.666 

85  X  lOp 

77.092 

130.909 

130.893 

132.296 

133.069 

90  X  lOp 

86.829 

147.284 

146.800 

148.152 

150.660 

95  X  lOp 

97.252 

164.881 

163.811 

164.694 

168.701 

100  X  lOp 

107.835 

183.740 

182.021 

181. 84o 

187.056 

105  X  105 

203.744 

201.498 

199.605 

206.266 

110  X  lOp 

223.790 

222.312 

217.973 

226.280 

115  X  10; 

244.778 

244.674 

237.744 

247.149 

120  X  lOp 

266.723 

268.3^ 

259.574 

268.807 

125  X  lOp 

289.673 

292.695 

282.934 

291.282 

130  X  lOp 

313.686 

316.602 

307.874 

314.606 

135  X  lOp 

338.867 

341.410 

334.455 

339.194 

ll*0  X  lOp 

365.240 

367.136 

362.782 

364.812 

ll*-5  X  lOp 

392.321 

393.809 

392.937 

391.695 

150  X  lOp 

419.986 

421.469 

422.555 

420.034 

155  X  lOp 

448.384 

450.164 

452.956 

449.700 

160  X  lOp 

477.191 

479.947 

483.937 

480.775 

165  X  10; 

506.834 

510.883 

513.  to8 

513.345 

170  X  lOp 

537.374 

542.993 

543.843 

547.574 

175  X  10; 

568.907 

575.236 

575.691 

580.072 

180  X  lOp 

606.519 

609.391 

607.563 

185  X  10; 

632.794 

639.119 

643.543 

642.936 

190  X  10; 

674.160 

676.304 

678.634 

195  X  105 

702.689 

709.215 

711.085 

712.460 

200  X  lOp 

738.494 

744.626 

747.457 

747.360 

205  X  lOp 

775.064 

783.186 

783.837 

210  X  lOp 

812.435 

819.664 

819.833 

822.144 

215  X  lOp 

850.412 

858.164 

857.787 

862.758 

220  X  lOp 

897.297 

906.208 

225  X  lOp 

930.670 

939.077 

939.348 

952.368 

230  X  10 

973.108 

981.827 

983.982 

998.529 
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Table  li)-  A  (cont’d) 


Velocity 

(ft/sec) 

0 

23.10 

42.16 

60 

80 

235  X  10^ 

240  X  lOp 
245  X  lOp 
250  xlOp 
255  X  10^ 

260  X  lOp 

265  X  lOp 
270  X  lOp 

275  X  lOp 
280  X  lOp 
285  X  lOp 

290  X  lOp 

295  X  lOp 

300  X  lOp 

305  X  lOp 

310  X  lOp 

315  X  lOp 

320  X  lOp 

325  X  lOp 

330  X  10; 

335  X  10; 

340  X  lOp 

345  X  105 

350  X  10^ 

355  X  10; 
360  X  lOp 
365  X  10^ 

370  X  10; 

375  X  10^ 
380  X  lOp 
385  X  lOp 
390  X  10; 
395  X  lOp 
400  X  10 

1017.023 

1062.885 

1107.591 

1152.523 

1197.995 

1244.175 

1292.193 
1342.545 
1394.287 
1447. 309 
1501.529 
1555.679 
1610.311 
1666.009 

1025.417 

1070.229 

1116.307 

1163.631 

1211.820 
1259.769 

1308. 649 

1358.935 

1410.496 

1463.583 

1518.384 

1574.651 

1629.161 

1684.821 

1029.857 
1077.403 
1126.683 
1174. 318 

1221.943 

1270.447 

1319.892 

1370-325 

1421. 859 
1474.498 
1528.308 
1583.374 
1639.790 
1697. 658 
1754.163 
1810.509 

1867.597 

1926.238 

1987.145 

2043.959 

2109. 716 
2171.228 

2233.716 

2297.411 

2362.016 

2428.009 

100 
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Table  Ito 

NORMAL  SHOCK  WAVE  CHARACTERISTICS  OF  A  TENTATIVE  VENUS  ATMOSPHERE 

TEMPERATURE  RATIO 


Velocity 

(ft/sec) 

Altitude  (Km) 

0 

25.10 

60 

80 

100 

- - j - 

20  X  lOg 

1.222 

1.1*07 

1.821* 

1.821* 

1.821* 

25  X  10^ 

1.585 

1.655 

2.270 

2.270 

2.270 

50  X  10^ 

1.566 

1.91*0 

2.776 

2.776 

2.776 

55  X  10^ 

1.776 

2.258 

5.51*0 

5.51*0 

5.51*0 

UO  X  lOp 

2.009 

2.618 

5.958 

5.958 

5.958 

h5  X  lOp 

2.261i 

5.012 

1*.61*6 

1*.61*5 

1*.61*5 

50  X  lot 

5.1*58 

5.585 

5.582 

5.570 

55  X  10^ 

5.895 

6.171* 

6.16I* 

6.156 

60  X  lOp 

1*.578 

6.997 

6.91*2 

6.862 

65  X  lOp 

1*.892 

7.81*7 

7.700 

7.505 

70  X  lOp 

5.1*00 

8.711 

8.1*58 

8.01*8 

75  X  10 

5.895 

9.551* 

9.208 

8.1*71* 

80  X  lOp 

6.567 

10.515 

9.977 

9.051* 

8.747 

85  X  10^ 

6.810 

11.051* 

10.716 

9.655 

8.796 

90  X  10^ 

7.208 

11.666 

11.1*02 

10.280 

8.753 

95  X  lOp 

7.551 

12.187 

12.006 

10.901 

8.931 

100  X  lOp 

7.989 

12.571 

12.1*91 

11.507 

9.368 

105  X  lOp 

12.866 

12.853 

12.087 

9.819 

110  X  lOp 

15.508 

15.066 

12.628 

10.277 

115  X  lOp 

11*.  157 

15.090 

15.122 

10.737 

120  X  lOp 

ll*.7l*2 

12.985 

13.516 

11.191 

125  X  lOp 

15.311 

13.015 

13.811* 

11.629 

150  X  lOp 

15.850 

13.1*97 

13.996 

12.045 

155  X  10 

16.281* 

13.980 

11*. 051* 

12.559 

lUO  X  lOp 

16.651 

11*.  1*57 

13.895 

12.628 

lh3  X  10 

17.092 

1U.923 

13.552 

12.851 

150  X  lOp 

17.61*1 

15.373 

13.521* 

15.019 

155  X  lOp 

18.251 

15.799 

13.532 

13.122 

160  X  lOp 

18.985 

16.195 

13.617 

15.152 

165  X  lOp 

19.71*0 

16.51*9 

11*.  317 

15.091 

170  X  lo| 

20.519 

16.867 

15.01*1* 

12. 908 

175  X  10 

21.517 

17.1*1*1 

15.792 

13.215 

180  X  lOp 

22.555 

18.1*92 

16.555 

14.749 

185  X  lOp 

25.522 

19.656 

17.516 

16.545 

190  X  10 

21*.  71*7 

20.951* 

18.856 

18.088 

195  X  lOp 

26.021 

22.291* 

20.1*1*0 

19.541 

200  X  lOp 

27.21*1* 

25.672 

21.936 

20.59I* 

205  X  lOp 

28.651 

25.095 

25.090 

21.224 

210  X  lOp 

50.229 

24.089 

21.804 

215  X  10 

51.556 

24.915 

22.104 

Table  14b  (cont'd) 

NORMAL  SHOCK  WAVE  CHARACTERISTICS  OF  A  TENTATIVE  VENUS  ATMOSPHERE 

TEMPERATURE  RATIO  Tg/T^ 


Velocity 


_  Altitude  ^Km 


23.10  i»2.l6  60 


220 

3( 

10: 

225 

X 

10; 

230 

X 

10: 

235 

X 

10: 

240 

X 

10; 

245 

X 

10; 

250 

X 

10; 

255 

X 

10: 

260 

X 

10^ 

265 

X 

10: 

270 

X 

10: 

275 

X 

10' 

280 

X 

10^ 

285 

X 

10^ 

290 

X 

10^ 

295 

X 

10^ 

300 

X 

10^ 

505 

X 

lOg 

310 

X 

10^ 

315 

X 

1O5 

320 

X 

10^ 

325 

X 

lof 

330 

X 

335 

X 

10§ 

340 

X 

345 

X 

102 

350 

X 

lOp 

355 

X 

lOg 

360 

X 

lOp 

365 

X 

lof 

370 

X 

lOp 

375 

X 

lOp 

380 

X 

10? 

CO 

X 

102 

390 

X 

lOp 

395 

X 

lOp 

400, 

X 

10^ 

32.660 

33.574 

34.313 

34.863 

35.216 

35.967 

36.811 

37.654 

38.495 

39.252 

39.774 

40.239 

40.656 

41.036 

41.585 

42.186 

42.770 


28.290 

29.002 

29.523 

30.085 

30.574 

30.990 

31.344 

31.707 

32.245 

32.727 

33.146 

33.493 

33.758 

33.928 

34.045 

34.638 

35.247 


25.543 

25.951 

26.108 

26.208 

26.187 

26.029 

26.319 

26.782 

27.230 


29.281 


40.126 


22.094 

21.734 

21.565 

21.942 

22.411 

22.883 


25.038 

25.220 

25.318 


26.206 


28.241 

29.117 


30.806 
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Table  14C 


N0RI4AL  SHOCK  WAVE  CHAHACTERISTICS  OF  A  TENTATIVE  VENUS  ATMOSPHERE 

DENSITY  RATIO,  p^/p^ 


Velocity 

Altitude 

('Km) 

(ft/sac ) 

0 

23.10 

42.16 

60 

80 

100 

20x10^ 

25x10^ 

30xl0; 

35x10^ 

40x10^ 

45x10^ 

50x10^ 

55x10^ 

60xl0^ 

65x10^ 

70x10^ 

75x10^ 

80x10^ 

85xio; 

90xio; 

95xlo; 

100x10^ 

110x10 

113x10^ 

120x1,0^ 

125xl0| 

130x10^ 

135x10^ 

140x10^ 

145x1 o; 
150xl0; 
i55xio; 
l60xl0^ 

165x10; 

170x10^ 

175x10; 

180x10^ 

185x10; 

190x10^ 

195x10; 

2ooxio; 

205x10; 

2ioxio; 

215x10; 

220xlo; 

225x10^ 

2.201 

3.089 

3.968 

4.772 

5.465 

6.134 

2.847 

3.844 

4.775 

5. 629 
6.353 
6.922 
7.490 
8.139 

8.  524 

9.064 

9.476 

9.792 

10.163 

10.654 

11.323 

12.137 

12.  370 

3.697 

4.724 

5.625 

6.421 

7.125 

7.675 

8.  097 
8.656 
9.095 
9.491 

10.  056 

10. 332 

10, 619 
10.984 
11.491 
12.191 
13.125 

14.212 

14.461 

14.681 

14.901 

15.154 

15.478 

15.914 

16.  512 
16.935 
17.057 
17.029 
16.742 

16. 443 
16.172 
15.970 
15.403 

14.943 
14.834 
14.893 
14.703 
14.505 
14.304 
14.053 

13.943 
13.961 

3.  697 
4.724 
5.625 
6.421 
7.125 
7.676 
8.103 
8.666 
9.171 
9.649 

10. 362 
10.680 
10.827 
10.932 

11.  067 

11.311 

11.728 

12. 340 
13.202 

14 , 509  . 

16,220 

17.651 

17.750 

17.804 

17.838 

17.879 

17.955 

18.093 

18. 318 
18.662 
19.131 
19.037 

17.971 

17 . 280 

17.319 

16.972 
16.477 
16.331 

16.130 

15.860 

15.696 

15.682 

3.697 
4.724 
5.625 
6.421 

7.125 

7.681 

8.122 

8.697 
9.272 
9.837 

10.  626 
11.692 
12.169 

12. 308 
12.253 
12.025 

11.655 
11.194 
10.695 

10. 522 

10.846 

11.393 

12.198 

13.337 

14 . 970 

17.377 

18.933 

20.418 

21.632 

20. 577 
19.800 
19.463 
19.726 
19.644 
18.408 

17.846 
17.660 
16.927 
16.302 
15.854 

15.656 
15.895 

11.568 

13.183 

15.171 

16.495 

16,840 

17.015 
17.032 
16,918 
16,703 
16.428 
16.139 
16.129 
16.223 
16.513 
17.089 
17. 903 
19.019 
20.556 
22.773 
22.658 
18. 651 
19.301 
19.590 
18. 466 
17. 622 
17.122 
17.025 
17. 394 
18.341 
20.010 
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Table  litc  (cont'd) 


NORMAL  SHOCK  WAVE  CHARACTERISTICS  OF  A  TENTATIVE  VENUS  ATMOSPHERE 

DENSITY  RATIO,  p^/p^ 


Velocity 

Altitude 

(Kk) 

(ft/sec) 

0 

23.10 

12.16 

60 

80 

100 

230x10^ 

11.083 

15.857 

16. 501 

21.100 

235x10; 

11. 308 

15.966 

17.191 

21.596 

2L0x10; 

11.695 

16.110 

18.096 

21.567 

2L5x10^ 

11.695 

16. 380 

19.258 

21.176 

250x10^ 

ll. 566 

16.677 

19.615 

21.375 

255x10^ 

11. 391 

16.910 

19.697 

21 . 277 

260x10“ 

11.203 

16.926 

I9.7I1 

21.200 

265x10^ 

11.175 

16. 9I6 

19.783 

21.167 

270x10^ 

11. 351 

17.027 

19.837 

21.202 

275x10; 

11. 581 

17.187 

19.931 

21. 3II 

280x10; 

11.818 

17.156 

20.062 

21.638 

285x10; 

15.II3 

17.871 

20.2I5 

22.171 

290x10; 

15.276 

18. 377 

2O.I99 

22.727 

295x10; 

15.3I8 

18.308 

20.816 

22,819 

300x10; 

15.110 

18.283 

21.  312 

23.093 

305x10; 

! 

21,156 

23. 319 

310x10; 

20.771 

23.311 

315x10; 

20. 39I 

22,776 

320xlo; 

20.180 

22,526 

325x10; 

20.216 

22.698 

330x10; 

20.289 

23.097 

335xl0; 

19.950 

22,805 

340xlo; 

19.623 

22.611 

3^5x10; 

19.332 

22.132 

350x10; 

19.108 

22.267 

355x10; 

18.891 

22.129 

360x10; 

18.761 

22.  O'?'! 

365x10; 

22.002 

370x10; 

21.951 

375xlo; 

21.867 

380xl0p 

21.795 

385x10; 

2I.7I7 

390x10; 

21.731 

395xlo; 

2I.7I9 

li-OOxlO'^ 

. 

21.761 
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Table  llvD 


NORMAL  SHOCK  WAVE  CHARACTERISTICS  OF  A  TEKTATIVE  VENUS  ATMOSEHERE 
VELOCITY  RATIO,  u^/u^ 


Velocity 

Altibuda 

(Km) 

(ft /sec) 

0 

23.10 

42.16 

60 

80 

20xlo| 

.454 

.351 

.270 

.270 

.270 

25x10^ 

.324 

.260 

.212 

.212 

.212 

30x10^ 

.252 

.210 

.178 

.178 

.178 

-J  P 

35x10^ 

.210 

.178 

.156 

.156 

.156 

14.0x10^ 

.183 

.157 

.140 

.140 

.140 

14.5x10^ 

.163 

.144 

.130 

.130 

.130 

50x10^ 

.134 

.124 

.123 

.123 

55x1o2 

.123 

.116 

.115 

.115 

6oxio^ 

.117 

.110 

.109 

.108 

65x10^ 

.110 

.105 

.104 

.102 

70x10^ 

.106 

.099 

.097 

.094 

75x1o2 

.102 

.097 

.094 

.  086 

80x10^ 

.098 

.094 

.092 

.082 

85x10^ 

.094 

.091 

.092 

.081 

.076 

90x10^ 

.088 

.087 

.090 

.082 

.  066 

95x10^ 

.082 

.082 

.088 

.083 

.061 

100x10^ 

.081 

.076 

.085 

.  086 

.059 

105x10^ 

.070 

.081 

.089 

.059 

110x10^ 

.069 

.076 

.093 

.059 

115x10^ 

.068 

.069 

.095 

.059 

120x10^ 

.067 

.062 

.092 

.060 

125x10^ 

.  066 

.057 

.088 

.061 

136x10^ 

.065 

.056 

.082 

.062 

135x102 

.063 

.056 

.075 

.062 

Il4-0xl0„ 

.061 

.056 

.067 

.062 

145x10^ 

.059 

.056 

.058 

.061 

150x10^ 

.059 

.056 

.053 

.059 

155x10^ 

.059 

.055 

.049 

.056 

l60xl0^ 

.060 

.055 

.046 

.053 

165x10^ 

.061 

.054 

.049 

.049 

170x10^ 

.062 

.052 

.050 

.044 

175x10^ 

.063 

.053 

.051 

.044 

l80xl0^ 

.065 

.056 

.051 

.054 

185x10^ 

.067 

.058 

.051 

.052 

190x10^ 

.067 

.058 

.054 

.051 

195x10^ 

.067 

.059 

.056 

.054 

200x10 

.068 

.061 

•  057 

.057 

205x10^ 

.069 

.061 

.059 

.058 

210x10^ 

.070 

.062 

.061 

.059 

215x10^ 

.071 

.063 

.063 

.058 

220x10^ 

.072 

.064 

.064 

.055 

225x10^ 

.072 

.064 

.063 

.050 

230x10^ 

.071 

.063 

.061 

.047 

235x10^ 

1 

.070 

1  .063 

.058 

.046 
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Table  Ito  (cont'd) 


NORMAL  SHOCK  WAVE  CHARACTERISTICS  OF  A  TENTATIVE  VENUS  ATMOSPHERE 
VELOCITY  RATIO,  u^/u^ 


89 


Table  I5 

EQUILIBRIUM  CONSTANTS,  (atm) 


T  (°K) 

Reaction 

1 

2 

-5 

5 

5 

xooo 

3.022"^® 

2.525"^^ 

1.46o4'5^ 

2.397'^® 

2.216“^^ 

1500 

1.215"^** 

1.375"^^ 

3.492"^^ 

5.054"^5 

5.438'^'^ 

2000 

8.970"33 

1.219"^^ 

6.332'^ 

2.664"33 

3.135"^ 

3000 

9.254"^ 

1.793 

1.855'^® 

1.9T9"^ 

2.519"^'^ 

4000 

2.828"^^ 

4.191"^^ 

2.108"^5 

2.765'^^ 

5000 

1.686"^ 

7.601"^° 

9.868“^^ 

3.290"^ 

6000 

6.088"^ 

1.209"^ 

2.995"^ 

3.984"’^ 

Table  I6 

EIDCIRON  CONCENTEIATIONS  FOR  PURE  CO„,  n  - 

2'  e 

(g-moles  electrons/^^.Oll  g  mixtiire) 
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Appendix  B 

THE  COMPUTATION  QF  TWQ-BIMEJISIOMAL  OBLIQUE- SHOCK- WAVE  CHARACTERTSTTCS 

FROM  NORMAT^  SHOCK- WAVE  DATA  — - 

Tbe  practical  problems  ustially  associated  -with  aerodynamics  generally 
involve  heat  transfer  and  aerodynamic  forces  vhich  arise  due  to  the  motion 
of  a  body  through  a  fluid.  In  the  present  application,  both  of  these  problems 
require  a  knowledge  of  the  so-called  invlscid  parameters  associated  with 
both  shock  waves  and  expansive-type  flows  such  as  the  Prandtl-Meyer  expansion 
around  a  comer.  The  purpose  of  this  appendix  is  to  furnish  a  practical  way 
to  calcuJ.ate  the  oblique-shock  parameters.  Since  aerodynamic  parameters 
change  only  in  the  direction  normal  to  shock,  only  the  oblique  shock  angle 

is  required  to  solve  the  problem,  ;Aien  the  normal-shock  solution  (Section  IV) 
is  knovm. 

The  present  method  may  also  prove  useful  for  counting  oblique- shock 
characteristics  for  any  selected  gas  composition  once  its  eqxiilibrlum  thenno- 
dynamic  properties  are  kiiown.  While  the  generally  used  iteration  scheme  for 
obtaining  normal-shock  data^^^  is  highly  convergent,  the  iteration  scheme 
for  oblique-shock  solutions  presented  in  Ref.  *20  is  not,  especially  near  the 
detachment  point.  It  is  presumed  that  if  it  is  practical  t©  obtain  normal- 
shock  data  for  a  particular  composition,  the  oblique- shock  characteristics 
can  be  obtained  by  this  method  for  a  particular  cace. 

As  previously  mentioned,  all  of  the  physics  of  the  problem  are  contained 
in  the  normal- shock-wave  data  and  the  fact  that  the  velocity  component  paral- 
lei  to  the  shock  wave  remains  constant  across  the  shock  (this  foUowa  directly 
frem  a  conservatlon-of-mass-and-momentum  analysis  across  a  flow  discontinuity). 
Iherefoie,  no  physical  assungrtions  about  the  flow  are  necessary. 
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What  is  usually  speelfled  In  an  oblique -shock  jaroblem  is  th<!  flow  de- 
fleetien  angle,  the  upstream  thersaodynamle  state,  and  the  upetream  velocity. 

If  the  -wave  inclination  angle  may  be  found,  all  downstream  parameters  may 
be  calculated.  Unfortunately,  no  explicit  expression  for  the  wave  inclination 
angle  exists,  and  a  practical  iteration  scheme  vdiich  always  converges  must 
be  found. 

Referring  to  Fig.  12,  the  following  trigonometric  relations  across  an 
oblique  shock  become  obvious; 


=  w^  sin  p 


(39) 


^2  _  tan  (p  -  9) 


tan  p 


(hO) 


The  same  conservation  equations  apply  to  the  normal  stream  cocixDnent  across 
a  flow  discontinuity  as  apply  to  the  stream  across  a  noirnal  shock  wave.  Since 
for  a  given  upstream  state  u^/u^  vs  u^  is  usually  plotted  for  a  particular 
gas  if  noimal  shock  wave  data  are  given,  this  relation  must  also  hold  for 
the  nomnal  stream  components  across  an  oblique  shock  if  the  upstream  theimo- 

^ame.  The  general  form  of  the  relation  is  shown  in 
quadrant  2  of  Fig.  I3.  qihis  curve  always  assumes  this  form  for  a  perfect 
gas.  Real-gas  effects  will  usually  destroy  the  monotonic  nature  of  this 
curve  and  will,  m  exceptional  cases,  cause  a  positive  slope  in  a  small  tu 

(20)  ^ 

range.  However,  this  presents  no  problem,  as  will  become  evident  later. 

Equation  (39)  is  very  simple  and  presents  no  cocrputatlonal  problem. 
Considering  Eq.  (kO),  the  following  equation  for  p  may  be  developed  after 


some  manipulation: 
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Fig,  13 — Graphical  solution  for  U2/U,  and  $ 
(weak-shock  solution) 


9^ 


Differentiating  Eq_.  (U-O),  the  maxirauiii  of  vs  P  with  0  a  fixed  paraioster 

oceurs  at 

P  =  +  I  (^2) 

and 


1  -  sin  0 
1  +  sin  0 


(1^3) 


Also,  referring  to  Fig.  I3,  it  is  seen  that  the  negative  sign  in  Eq.  (4l) 
corresponds  to  0  <  p  <  6/2),  and  the  positive  sign  corresponds  to 

(45°  +  0/2)  <  P  <  90°*  Note  that  0°  <  p  <  90°  is  the  only  range  of  interest. 
For  reasons  that  will  hecoaae  apparent  later,  it  is  usually  p  ^Ich  raust  he 
calculated  from  a  known  u^/u^,  rather  than  vice  versa.  Since  Eq.  (4l)  re¬ 
quires  a  rather  tedious  cosnputatlen  if  done  hy  hand  in  an  iteration  sclieme. 
Fig.  Ik  is  Included  for  0  <  0  <  30°. 

It  may  be  shown  that  p  <  (45°  +  0/2),  i.e.,  P  <  P  corresponding  to 

(u  /u_  )  for  a  given  0,  must  always  be  in  the  weak- shock- wave  region, 
majc 

Stating  this  another  way:  if,  for  a  given  \ipstream  state,  P  were  plotted 
vs  0  with  upstream  velocity  as  a  parameter,  p  <  (4-5°  +  6/2)  would  always 
be  less  than  the  p  dictated  by  the  line  connecting  the  mp-rimum  0  solutions 
for  a  given  upstream  velocity.  Ihls  may  be  seen  for  a  perfect  gas  by  i^ottlng 
a  line  p  ®  (45°  +  Q/S)  on  pp.  42  and  43  of  Ref.  22.  The  intersection  of  this 
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Fig. 14 — Graphical  solution  for  shock-wave  angle 


line  with  0  =■  0  at  Mach  =  »  Is  not  a  coincidence,  hut,  as  may  be  shown, 
max 

a  consistent  result.  Uae  existence  of  real-gas  effects  does  not  alter  the 
original  statement,  and  as  the  upstream  velocity  becomes  large  p  =  (45°  +  6/2) 
approaches  the  limit  of  the  weak- shock  solutions.  At  low  upstream  Mach 
nttmbers  a  large  range  of  weak-shock  solutions  lies  where  p  >  (45°  +  9/2),  as 
may  be  seen  on  p.  42  of  Ref.  22  (at  low  Mach  numbers  a  gas  usually  behaves 
as  a  perfect  gas).  At  any  rate,  the  najority  of  physically  possible  solutions 
where  the  upstream  Mach  number  is  greater  than  2  will  occur  where  P  <  (45° 

+  6/2) j  this  will  turn  out  to  be  fortunate,  considering  convergence  of  the 
method. 

Now  a  solution  satisfying  Eqs.  (39)  and  (40)  and  the  normal  shock 
relation  of  xig/tL,  vs  for  a  given  upstream  state,  w^  and  0,  is  desired. 

A  qualitative  plot  is  now  constructed  in  Fig.  I3  as  follows: 

<cjuadrant  1:  A  functional  relationship  between  and  p  for 

the  given  0  (Eq. (40)) 

Quadrant  2:  The  normal-shock  relation  u^/u^  vs  u^  for  the  given 
upstream  state  (contains  the  physics  of  the  problem) 

Quadrant  3:  A  functional  relation  between  u^  and  p  for  the  given 

^1  ^  ^1 

(^adrant  4:  A  45°  line  to  graphically  convert  the  p  axis  in 
quadrant  1  to  the  P  axis  in  quadrant  3 
Line  ABCDA  can  be  seen  to  represent  the  wak-shock  solution.  The  procedure 
to  arrive  at  this  solution  is  as  follows: 

1.  Assume  p  =  90°  in  quadrant  3  such  that  ^  =  '^2. 

2.  From  u^  =  w^  find  Ug/iij,  quadrant  2  (point  F). 

3.  From  Ug/u^,  find  p  in  quadrant  1  (point  G). 
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4.  Carry  over  this  p  to  quadrant  3  throu^  quadrant  14-  (point  H).  It 

is  now  apparent  that  throu^  the  monotonic  nature  of  all  the  functions  the 

solution  has  been  bracketed  between  p  =  90°  and  p_  = 

G  H 

5-  Obtain  u^  in  quadrant  3  (point  J). 

6.  Obtain  in  quadrant  2  (point  K). 

7.  Obtain  3  in  quadrant  1  (point  L).  Again  it  is  obvious  that  the 

solution  is  now  bracketed  between  and 

Li  G  H 

8.  Continuing  in  this  clockwise  fashion,  it  is  seen  that  the  process 
is  rapidly  convergent. 

It  is  important  to  note  that  any  attempt  to  construct  this  solution 

in  a  coiinterclockwise  manner  would  have  failed,  since  the  process  would 

diverge.  It  is  easily  seen  that  this  process  converges  In  the  same  manner 

■vdien  the  normal- shock  relation  in  quadrant  2  contains  inflection  points 

but  retains  negative  slope.  It  will  be  stated  without  demonstration  that 

this  process  also  converges  when  the  normal-shock  relation  contains  a  small 

region  of  positive  a]x)pe  as  previously  mentioned.  In  this  case,  however, 

the  convergence  to  the  solution  is  from  one  side  only;  the  solution  is  not 

continually  being  bracketed  as  above. 

One  other  difficulty  may  be  encountered.  On  the  second  iteration 

point  K  may  correspond  to  a  Ug/u^  >  (ug/u^)  .  In  this  case  the  starting 

max 

point  for  the  third  trial  is  the  p  corresponding  t®  (Ug/iL,  )  ,  l.e., 

p  XI  (h5  +  9/2).  If  on  the  third  trial  the  ratio  Ug/u^  obtained  la  quadrant 

2  Is  greater  than  (vig/u^)  ,  the  solution  must  be  \diere  p  >  (1^5°  +  0/2). 

max 

Figure  15  Illustrates  the  case  mentioned  immediately  above.  Following 
the  lines  JKLMK  and  OPQR,  it  is  appeor^nt  that  the  solution  must  not  be 
where  P  <  (h5°  +  6/2).  Now  starting  again  at  point  J  and  proceeding  to  K, 
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continue  the  KL  line  to  S  ^ere  3  >  (45°  +  e/2).  Continue  with  a  clockwise 

iteration  starting  from  JTCSTV,  always  selecting  p  >  (45°  +  0/2)  in  qttadrant 

1;  -this  is  seen  to  converge  to  the  solution  ABCDA.  Alternatively,  starting 

from  (ug/uj^)  at  point  0  and  proceeding  counterclockwise .  the  iteration 
max 

begun  by  OWXYZ  is  seen  to  be  converging  to  the  solution  EFGHE,  It  may  also 
be  seen  that  if  p  is  originally  selected  such  tlsat  p„  <  p  <  p  ,  the  clock- 
ivise  iteration  will  converge  to  ABCDA  and  coimter clockwise  iteration  will 


converge  to  EFGHE.  It  will  also  be  stated  without  demonstration  that  these 
conclusions  are  not  altered  if  small  regions  of  positive  slope  occur  in  the 
normal-shock  relation  (quadrant  2).  It  then  becomes  apparent  that  EFGHE  is 
the  weak-shock  solution  and  ABCDA  is  the  strong-shock  solution. 

With  all  cases  now  considered,  the  general  method  to  be  followed  in 
obtaining  a  solution  is  outlined  below. 

Given:  Upstream  state,  0,  w^^,  and  the  normal- shock  relation  between 
and  u^  for  the  gi-'/en  upstream  thermodynamic  state. 

A.  Weak- shock  solution 

1.  Assucme  p  =  90°,  i.e.,  u^  =  w^ 

2.  From  the  normal  shock  relation  obtain  u^/u^ 

3.  From  Eq.  (4l)  or  the  included  graphs,  obtain  a  new  p,  choosing 

P  <  (45°  +  0/2) 

4.  From  this  p  compute  a  new  u^  from  Eq.  (39) 

5.  From  the  new  u^  and  the  normal- shock  relation  obtain  a 
new  Ug/u^ 

6.  If  this  Ug/u^  <  (ug/u^)  for  the  given  0  from  Eq.  (40), 

max 

repeat  this  procedure  to  the  desired  degree  of  convergence 

7.  If  this  ^  assume  p  =  45°  +  0/2  and  connjute  a 

max 
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new  fram  Eq,.  (39) 

8,  Using  "this  and  the  Kormal-shoek  relation  obtsiln  a  new 

9.  If  this  <  (ug/u^)  proceed  as  in  step  6 

Kiax 

10.  If  this  Ug/xij^  >  (u^/ii^)  ,  the  solxition  will  not  lie  xdxere 

o 

3  <  (45  +  0/2);  assxme  u^/xx^  =  (ug/xj^)  ;  obtain  xu  from 

max 

the  normal  shoek  relation 

11.  Compute  3  from  Eq.  (39)  and  this  xj^ 

12.  Using  tills  p  compute  u-^/u^  from  Eq.  (40)  or  obtain  p  from 
the  included  charts 

13-  Using  this  obtain  u^  from  the  normal-shock  charts 

l4.  Repeat  steps  11  —  I3  until  the  desired  degree  of  convergence 
is  reached 

B.  Strong-shock  solution 

1.  Assume  p  =  90°,  i.e.,  ^  '^1 

2.  From  the  normal  shock  relation  obtain 

3.  From  Eq,  (4l)  or  the  Included  charts  obtain  a  new  p,  choosing 
the  p  >  (45°  +  0/2) 

4.  Compute  a  new  u^  from  Eq.  (39) 

5.  Repeat  steps  2-4  until  the  desired  degree  of  convergence 
is  reached 

It  may  be  shox/n  that  if  the  normal  shock  solution  is  obtained  by  an 
iteration  scheme  that  Is  stopped  when  two  consecutive  valxies  of  u^/u^  differ 
by  no  more  than  G.OQOl,  a  consistent  degree  of  convergence  in  this  method 
is  arrived  at  when  txro  consecutive  values  of  p  differ  by  no  more  than  O.O5 
per  cent,  i.e.,  2sp/p  <  0.0005;  however,  this  is  by  no  means  absolutely  true 
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over  the  entire  range.  The  rapid  convergence  of  this  method  is  illustrated 
by  three  numerical  examples  for  air  given  below. 

It  should  be  kept  in  mind  that  this  procedure  is  not  particularly 
useful  uhen  a  cor '.pie  te  set  of  oblique -shock  charts  for  a  kno'vm  atmosphere 
is  to  be  constructed.  In  such  a  case  1  that  needs  be  done  is  to  specify 
a  p  and  iterate  the  normal-shock  equations.  The  corresponding  9  can  -then 
be  immediately  calculated. 


EXAMPLES 


1.  Very  Weak  Shock 


Gi-ven: 


w^  =  8000  ft/sec  in  air 

Altitude  =  82,3^5  ft  (1956  ARDC  model  atmosphere ) ^ 


9  =  20^ 


Find:  p  for  the  -weak- shock  solution 

Solution:  (a)  Assvune  p  =  90°/  then  u^  =  8OOO  ft/sec 

From  the  normal-shock  data  in  Ref.  20,  u^/u^  =  0.1^35 
From  Eq.  (4l),  taking  the  negative  sign  for  p  <  (45°  +  e/2) 
P  =  23.58° 

The  solution  is  now  bracketed  between  23.58°  and  90° 

(b)  Taking  p  =  23*58°  and  using  Eq.  (39)/  u^  ■=  3201  ft/sec 
From  the  normal- shock  data  =  0,2407 

From  Eq.  (4l),  p  =  26.98° 

The  solution  is  now  bracketed  between  26.98°  and  26.16° 
Contln'uing: 

(b.l)  =  36^)  ft/sec 


0.2198 
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p  =  26.16° 

(h.2)  =  3527  ft/  sec 

=  0.2243 
P  =  26.33° 

(b.3)  =  3548  ft/  sec 


—  =  0.2230 
P  =  26.29° 

(b.4)  =  3544  ft/  sec 


— -  =  0.2230  within  the  accxiracy  with  lAiich  the  charts 
can  be  read 

P  =  26.29°  vdiich  is  the  solution  compared  with 
P  =  26.3°  obtained  from  Ref.  20 

— — g.^yoi^-Sh°ck  Solution  and  Illustration  of  Convergence  Difficulties  for 

Gi-ren  =  3000  ft/sec  in  air 

Altitude  =  175, 348  ft 

9  =  30° 

Find:  p  for  the  strong-shock  solution 

Solution:  First,  hunting  for  the  weak-shock  solution,  assume 

P  =  90° 


=  3000  ft/sec 

From  the  normal-shock  data 

=■  0.2740 
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From  Fig.  14  for  9  =  30°, 


0.2740,  and  p  <  (45°  +  0/2) 


P  =  45.67° 

From  Eq.  ( 39 ) 

=  2176  ft/sec 

and  from  the  normal -shock  chart 


u. 

X 


=  0.3866 


But,  as  can  he  seen  in  Fig.  l4,  this  is  greater  than 

(?) 

\  /max 

Therefore,  assume 

\^L 

From  Eq.  ( 39) 

u^  =  2598  ft/sec 


P  =  (45°  +  ■^)=  60° 


and  from  the  shock  chart 
^2 

■—  =  0.3140 

’^2 

vhich  is  now  less  than  —  and  the  process  will  con- 

max 

verge  to  the  weak-shock  solution  vhere  p  <  (45°  +  9/2). 
Looking  for  the  strong-shock  solution,  choose  p  =  90°, 

=  3000  ft/sec.  Then  as  before  u^/u^^  =  0.2740.  Now 
from  Fig.  l4  pick  p  >  (45°  +  0/2), 

P  =  74.32° 

From  Eq.  (39) 

=  2888  ft/sec 
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Fr®ni  "the  ahoek  cfasurfcs 
Up 

=  0.2843 
From  Fig.  14 

P  =  73.07° 

Contliming: 

u,  =■  2870  ft /sec 
Up 

-f  =  0.28i^8 
P  =■  72.99° 

=  2869  ft/sec 

TT  =■  0.2850 
P  =  72.96^ 

And  since  en  the  last  trial 


0.00041  <  0.0005 


we  step  the  process  and  taJie  the  strong ^hock  solution  as 
P  =■  72.96° 


3.  Weak-aiock  Solutian  on  Right  Half  of  u^/u^  vs  p  Curve 

Given:  Mach  number  =  I.9  in  air  which  is  assumed  a  perfect  gas 

Noirmel-sheck  relations  are  taken  from  Ref.  22,  0  =  20° 

Find:  p  for  weak  shock 

Solution:  Clearly  assuming  p  =  90°,  =  1.9  vill  yield 

normal 
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Contintiing: 


=  1.608 

*^n 

p  =  57.81° 
~  =  0.488U 
=  1.6091 

n 

p  =  57.87° 
—  =  0.4883 
=  1.6094 
p  =  57.89° 

and  since 


^  ^  57.89  -  $7.87 
p  57r^9 


0.0003  < 


the  process  is  stopped,  and 
P  =  57.89° 


0.0005 


compared  to  the  solution  in  Ref,  22  of 
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Appeadix  C 

CALCXJLAnon  OF  ELiECJTRON  CONCaSMTEATIONS  AT  LOW  TEMPERA'PURES  TK  THTB 

C-M-0  GASEOUS  SYgTEM  “ 

It  is  kncnm  that  strong  interference  vlth  .radio  or  telenetiy  trana- 
mission  from  ataaospheric-entry  vehicles  is  possible  if  the  transmiasion 
frecLuen<gr  is  equal  to  or  less  tham  the  plasma  frequent^  of  the  sur¬ 

rounding  the  antenna  and  if  the  collision  frequency  is  large  enou^.  Plasma 
freqiKnoy  is  directly  dependent  upon  the  square  root  of  the  electron  con¬ 
centration  in  the  surrounding  medium,  and  it  may  be  desired  to  know  this 

-17 

concentration  to  as  low  as  10  g-mole  of  electrons  per  g-mole  of  gas 
mixttire.  Kiis  concentration  corresponds  to  a  transmisBion  frequency  of  10 
Me  throu^  a  medixim,  for  examjjle,  at  a  pressure  of  10^  atm  and  a  tenqperature 
of  1000°K  in  the  present  tentative  Venus  atmasihere. 

Since  the  equilibrium  con^wsltion  of  gases  coii5Juted  by  the  method  of 
Ref.  18  does  not  show  the  trace  constituents,  it  is  not  suited  to  the  cal¬ 
culation  of  electron  concentrations  at  low  ten^^ratures  if  the  major  con¬ 
stituents  of  the  gas  mixture  are  known.  In  particular,  the  present  calcu¬ 
lations  apply  to  the  C-N-0  gaseous  mixture  but  the  general  method  may  be 
applied  to  any  arbitrary  mixture.  Hie  calculated  electron  concentration  is 
based  on  an  ideal  gas  mixture  in  thermodynamic  equilibrium  when  the  mjor 
constituents  of  the  gas  mixture  are  known.  Since  the  major  constituents  of 
the  gas  are  generally  computed  on  the  basis  of  an  ideal  gas  in  thermodynamic 
equilibrium,  these  conditions  are  consistent. 

For  a  chemical  reaction  of  ideal  gases  in  thermodynamic  eqioilibrium, 
the  equilibrium,  constant  may  be  inrltten 
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•i*ere  1  denotes  the  products,  J  the  reactants,  emd  A  the  change  as  measured 
by  products  minus  reactants,  an  examination  of  the  reaction  energies 

for  the  C-N-0  system  in  Table  4,  it  is  evident  that  at  low  ten5>er  at  vires 
the  electrons  may  to  a  close  approximation  be  considered  to  be  produced 
solely  by  the  following  reactions: 

1.  COg  ^  COg'*'  +  e"  AE°  =  317,731  cal/g-mole^^^^ 

2.  CO  ^  CO  +  e  AE°  =  323,180  cal/g-mole^^^ 

3-  Og  —  Og'*'  +  e"  AE°  =  277,900  cal/g~mole^^^ 

4.  0  ^  0  +  e  AE°  =  314,048  cal/g-aole^^^ 

5*  C  c  +  e  AE°  =  259,842  cal/g-mole^^ ^ 

6.  NO  KO"^  +  e"  AE°  =  213,400  cal/g-mole  ^  ^  ^ 

Using  the  le.-w  of  partial  pressures  the  equilibrium  constant  for  the  corre¬ 
sponding  reactions  may  be  expressed  as 
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(47) 

(48) 

(49) 

(50) 


Since  it  is  assumed  that  all  of  the  electrons  cone  from  the  above 
ionized  species 


^  ”cOK  +  +  +  nj^O+  (51) 

Combining  Eqs.  (t5)  -  (51),  the  electron  concentration  on  the  basis  of  n 
g-moles  of  the  mixture  is 


Some  equilibrium  constants  have  been  calculated  and  are  shown  in  Table  I5  for  a 
range  of  1000°K  to  6000°K.  The  amounts  n^  of  the  major  constituents  are  assumed 
known  and  may  be  placed  In  the  conq>utational  Eq.  (52)  along  with  the  equili¬ 
brium  constants,  and  the  electron  concentrations  may  be  determined.  It  is 
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suggested,  howsver,  that  the  accuracy  obtained  by  this  method  be  limited  to 
two  significant  figures  at  temperatures  where  electron  concentration  be¬ 
comes  significant,  since  the  molar  constituents  assumed  in  the  calculation 
from  Eq.  (52)  do  not  account  for  the  presence  of  electrons. 

For  Illustrative  purposes,  results  of  some  of  these  computations  for 
pure  COg  corresponding  to  those  of  Ref.  ( I7 )  have  been  made  and  are  pre¬ 
sented  in  Table  I6. 
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